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Abstract 

Four-dimensional asymptotically flat spacetimes at spatial infinity are denned 
from first principles without imposing parity conditions or restrictions on the Weyl 
tensor. The Einstein-Hilbert action is shown to be a correct variational principle 
when it is supplemented by an anomalous counter-term which breaks asymptotic 
translation, supertranslation and logarithmic translation invariance. Poincare trans- 
formations as well as supertranslations and logarithmic translations are associated 
with finite and conserved charges which represent the asymptotic symmetry group. 
Lorentz charges as well as logarithmic translations transform anomalously under 
a change of regulator. Lorentz charges are generally non-linear functionals of the 
asymptotic fields but reduce to well-known linear expressions when parity conditions 
hold. We also define a covariant phase space of asymptotically flat spacetimes with 
parity conditions but without restrictions on the Weyl tensor. In this phase space, 
the anomaly plays classically no dynamical role. Supertranslations are pure gauge 
and the asymptotic symmetry group is the expected Poincare group. 
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1 Two puzzles in asymptotically flat spacetimes 



Both in Hamiltonian framework or in covariant phase space methods, asymptotically flat space- 
times at spatial infinity have only been defined when parity conditions on the first order part 
of the boundary fields are imposed. These conditions have been introduced in [1] in the Hamil- 
tonian formalism and later in the covariant phase space, see [2], so that Lorentz charges are 
finite and so that the canonical structure, or in Lagrangian formalism the covariant symplectic 
structure, is also finite. Now, it was noticed in [3] (see also earlier work [4, 5, 6, 7, 8, 9, 10, 11]) 
that asymptotically flat spacetimes at spatial infinity admit a variational principle whether or 
not parity conditions on the first order part of the boundary fields in the asymptotic cylindrical 
or hyperbolic radial expansion hold, at least when one neglects boundary terms at the past and 
future boundaries. 

It is then intriguing that even though the action is finite, the symplectic structure and 
conserved charges are infinite when parity conditions do not hold. One might think that the 
action determines the entire dynamics and therefore in particular the symplectic structure and 
conserved charges but this expectation does not seem to be realized. This constitutes the first 
puzzle that we will resolve in this work 1 . 

Asymptotically flat spacetimes are defined by boundary conditions at infinity. The class 
of diffeomorphisms which preserve the boundary conditions are the allowed diffeomorphisms. 
Allowed diffeomorphisms associated with non-trivial conserved charges - the large diffeomor- 
phisms - modulo diffeomorphisms associated with zero charges - the pure gauge transformations 
- define the asymptotic symmetry group. For asymptotically flat spacetimes at spatial infinity 
when parity conditions are imposed, the asymptotic symmetry group has been worked out both 
with Hamiltonian and covariant phase space methods. In Hamiltonian framework, it is just the 
Poincare group. In that framework, parity odd supertranslations also preserve the boundary 
conditions but they are associated with vanishing Hamiltonian generators [1]. In covariant phase 
space, the asymptotic symmetry group has been shown to be the Poincare group only in the 
truncated phase space where part of the first order fields have been set to zero, by setting the first 
order magnetic part of the Weyl tensor to zero. This last condition also fixes all supertranslations 
[5, 15]. The equivalence of Hamiltonian and Lagrangian formalisms tells us that there should 
exist a covariant phase space which allows for more generic first order fields (see also comments 
in [16]) and such that the covariant equivalent of the Hamiltonian odd-supertranslations act as 
pure gauge transformations. We will show that indeed such a generalization of the covariant 
phase space with parity conditions exists. 

The existence of such a consistent covariant phase space or canonical phase space does 
not shed a complete light on the choice of boundary conditions. In particular, it is not clear 
whether or not covariant phase spaces related to each other by unallowed diffeomorphisms are 
inequivalent. This constitutes our second puzzle. Both in the usual covariant phase space and 
in Hamiltonian framework, different logarithmic translations or generic parity supertranslations 
are not transformations associated with finite Hamiltonian or Lagrangian generators when one 

1 We thank D. Marolf and A. Virmani for drawing our attention to this issue and A. Ashtekar for emphasizing 
the role of past and future boundary terms in the variational principle. A hint that parity conditions might not 
be required in a phase space including counter-terms was provided in [12] where it was shown that counter-term 
charges [3, 13, 14] and Ashtekar-Hansen charges [5] are equivalent without using parity conditions, apart for the 
four lowest harmonics of the first order field to fix logarithmic translations, in a restricted phase space where 
supertranslations are also completely fixed. 
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uses the standard canonical bracket or symplectic structure. There is therefore no covariant 
phase space or canonical space that encompass spacetimes or initial data surfaces with such 
different asymptotic transformations. Now, either these transformations are unphysical or they 
are physical. On the one hand, if they are unphysical, one would expect that they are pure gauge 
(see [17] for arguments that logarithmic translations are unphysical). The fact that logarithmic 
translations and generic parity supertranslations are not degenerate directions of the symplectic 
structure - they are not allowed directions of the symplectic structure in the first place - is 
however in tension with the intuition that pure gauge transformations should be degenerate 
directions of the symplectic structure. On the other hand, if these transformations have a 
physical content, there should be a way to regularize the infinities. If an enlarged phase space 
exists where both logarithmic and generic supertranslations are allowed in the first place, it would 
allow to settle these questions. However, there does not exist in the literature a construction of 
a consistent phase space where parity conditions have not been imposed. 

These two puzzles have in common that they rely on the standard covariant phase space 
symplectic structure or in Hamiltonian formalism on the canonical bracket defined from the 
canonical fields 3 gij and yr ij ". Now, it is important to remember that the algebraic derivation of 
the covariant phase space symplectic structure from the bulk Lagrangian suffers from ambiguities 
[18, 19]. Also, the definition of the canonical structure depends on what fields are considered to 
be canonical. For example, in [1], it was pointed out that the asymptotic values for the shift and 
lapse functions should be considered as additional canonical variables in addition to s gij and 
7r* J . The consideration of additional canonical variables might then lead to a modification of the 
canonical structure. It is also important to observe that in the arguments of [3] the variational 
principle is defined only under the assumption that future and past boundary terms can be dealt 
with without affecting the analysis at spatial infinity. 

2 Main results and outline 

The key technical result in this paper is the fixation of ambiguities in the symplectic structure 
and the canonical bracket. We show that the ambiguities can be fixed from first principles by 
requiring the existence of a variational principle taking into account boundary terms at past and 
future infinity in addition to spatial infinity. Our definition of symplectic structure then amounts 
to the prescription given in [20] which was developed following the counter-term methods [21, 22, 
23, 24, 25, 10] in the framework of the AdS-CFT correspondence [26, 27] . An important difference 
with earlier treatments of asymptotically flat spacetimes is the inclusion of a translation breaking 
counter-term, a translation anomaly, in the action. As a consequence of the addition of boundary 
counter-terms to the canonical structure and symplectic structure, no parity conditions are 
required at any step in the construction of the phase space and finite conserved charges can 
be defined from first principles. As a result, supertranslations and logarithmic translations 
are associated with non-trivial conserved charges. These conserved charges are trivial when 
parity conditions hold. Our Poincare conserved charges reduce to the standard ADM [28, 1] or 
Lagrangian charges [29, 5, 30, 31, 32, 18, 19, 33] when parity conditions hold. We also define a 
covariant phase space of asymptotically flat spacetimes with parity conditions that generalizes 
previous constructions [2, 34]. In this phase space, the anomaly does not contribute to the 
symplectic structure so most of the conventional properties of asymptotically flat spacetimes are 
kept, e.g. the Poincare group is the asymptotic symmetry group. These and other consequences 
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of our construction will be addressed in the following sections. 

We start in Section 3 by providing the detailed boundary conditions we will use throughout 
the main text. We make connection with previous choices of boundary conditions and indicate 
how they are generalized. In Section 4, we provide a first-principle derivation of the asymp- 
totic dynamics of Einstein gravity, with our flat boundary conditions, starting from the action 
principle and deriving the asymptotic equations of motion, the symplectic structure, the co- 
variant phase space charges and the asymptotic symmetry group. In Section 5 we present the 
main milestones on how our construction fits in the Hamiltonian framework. We end with a 
summary and detailed discussion in Section 6. Some additional details on the boundary con- 
ditions and the equations of motion are relegated to appendices. We provide in Appendix A 
a comparison between covariant and 3+1 boundary conditions. In Appendix B we obtain a 
classification of specific tensor fields useful to characterize the algebra of tensors at second order 
in the asymptotic expansion when all first order fields are considered. In Appendix C we derive 
the asymptotic equations of motion in the hyperbolic representation of spatial infinity up to 
second order in the radial expansion. 

We refer the reader to the companion paper [12] for a review of the hyperbolic representation 
of spatial infinity, our notations, and the presentation and proofs of some lemmas and properties 
of tensor fields on the hyperboloid. Although this paper deals only with smooth fields, a supple- 
mentary lemma on singular tensor fields on the hyperboloid relevant to describe NUT charges 
is stated in Appendix D for completeness. 

3 Specification of boundary conditions 

The boundary conditions are part of the specification of the theory. If the boundary conditions 
are too strong, there is no interesting physics in the restricted phase space. If they are too loose, 
no phase space can be defined at all because physical quantities like the total energy or the 
angular momenta diverge. Starting from a set of boundary conditions, one can try to enlarge 
them such that all physically reasonable solutions of the equations of motion lie inside the phase 
space and such that all physically relevant quantities like action, symplectic form and conserved 
charges are well-defined. Note that there is no theorem guaranteeing that this procedure leads 
to a unique "largest" phase space. In the sequel, we present boundary conditions in covariant 
phase space and 3+1 formalism which are more general than the ones discussed for example 
in [1] and [2]. The consistency of these conditions will be discussed in the later sections. The 
relationship between boundary conditions in the canonical and covariant formalisms is further 
discussed in Appendix A following earlier comparisons [15, 16]. 

3.1 Boundary conditions in the hyperbolic representation of spatial infinity 

In the hyperbolic representation of spatial infinity proposed in [5] and developed in [35, 36, 32], 
one demands that there should exist a coordinate system (p, r, 6, 4>) such that the four-metric 
takes the form 

ds 2 = (l + — + otp-^dp 2 + 0(p°) dpdx a + (p 2 /^ + ph$ + o(p l )y x a dx b , (3.1) 

where n£j} is the unit hyperboloid 

ds^ = hf l jdx a dx b = -dr 2 + cosh 2 r(d6 2 + sin 2 Ode/) 2 ), (3.2) 
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and a and h ab can be considered as a scalar and a tensor field defined on the hyperboloid. 
Minkowski spacetime takes this form when one transforms the standard spherical coordinates 
into p 2 = rjp LV x^x v = r 2 — t 2 and tanhr = t/r. We take as our covariant phase space all metrics 
that can be put into the form (3.1) and where k ab , defined as 

k ab = hV + 2oh ( $, (3.3) 

obeys the constraints given below in (3.6). In the following, D a denotes the covariant derivative 
associated with the metric on the hyperboloid h^? and □ = V a T> a . Indices are raised with 
h^ ab . 

As pointed out in [35, 36], the phase space is invariant under Poincare transformations, 
supertranslations and logarithmic translations. Logarithmic translations are defined as p — > 
p + log pH(x a ) + o(p°), x a -> x a + p- l {\ogp + 1) V a H{x b ) + o{p- 1 ) where H ab + h$H = 2 . 
They act on the first order fields a and k a b as 

5 H a = H, 5 H k ab = 0. (3.4) 

The transformations with H arbitrary are the so-called logarithmic supertranslations that could 
be used to set a = at the expense of introducing a logarithmic term of the form plog p(H ab + 
h^jj H)dx a dx b in the metric. We choose to set this logarithmic term to zero in (3.1) in order 
to follow the notation of previous literature. Logarithmic supertranslations are therefore not 
allowed diffeomorphisms with our choice of boundary conditions 3 . 

Supertranslations are defined as the transformations p — > p + uj(x a ) + o(p°), x a — > x a + 
p— 1 X> a oj(x b ) + o(p~ 1 ). When acting on the leading order fields a and k ab , we find 

6 U (7 = 0, 5 w k ab = 2(V a V b co + uh!$). (3.5) 

The trace of k ab , k = h^ ab k ab is unconstrained by the equations of motion derived from 
Einstein's equations while the divergence of k ab has to obey T> b k ab = T> a k. Supertranslations 
that change the trace k are given by transformations (3.5) that obey (□ + 3)cj ^ 0. As we will 
show in Section 4.5, one can associate charges to such transformations that are not conserved at 
spatial infinity. Since we aim at defining a phase space with only conserved charges at spatial 
infinity, we will discard such transformations by fixing the trace of k to an arbitrary scalar k. 
Consistently with the equations of motion, we impose the boundary conditions that k ab has fixed 
trace and divergence 

k = fc(°> ab k ab = k, V b k ab = V a k. (3.6) 

For simplicity, we will fix k = 0. As we will see in Section 4.1, the conditions (3.6) turn out to 
be sufficient in order that the Mann-Marolf action [3] be a variational principle on the phase 
space. Accordingly to the restriction (3.6), we restrict the function u;(x a ) on the hyperboloid to 
obey 

(□ + 3)cj = 0, (3.7) 

2 Our definition is equivalent to the linear combination of a logarithmic translation as defined in [35, 36] with 
a translation. 

3 Note that our second puzzle also applies to logarithmic supertranslations but for simplicity we will not consider 
these transformations in this work. 
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which is a severe restriction with respect to the definition of supertranslations used e.g. in [2]. 
Translations are singled-out as the four transformations obeying 

u a b + hf b 'u = 0, Lo ab = V a V h uj . (3.8) 

By convention, we will refer to supertranslations as the transformations generated by oj which 
obey (3.7) but do not obey (3.8). 

In order to make connection with previous treatments, let us review what the boundary 
conditions of [5] impose in addition to (3.1). As we review in Appendix C, once fixing k = 0, 
Einstein's equations imply that 

(curl (curl kj) = (□ - 3)k ab = 0, (curlT) ab = e a cd V c T db . (3.9) 

V / ab 

It has been shown in [5] (referred to as Lemma 1 in [12]) that these equations imply that k ab is 
determined in terms of a scalar /3 as 

cur\k ab = V a V b (3 + h { ° b ) p , (3.10) 

with (□ + 3)/3 = 0. In [5], it is imposed that 

/3 = 0. (3.11) 

Once the curl of k ab is zero, one can use again the Lemma in [5] to express k ab = V a V b I + h$I 
with (□ + 3)/ = 0. One can then fix completely the supertranslations by imposing 

k ab = 0. (3.12) 

In addition to this, it can be imposed (see e.g. [2] even thought this condition was not 
imposed in the original reference [5]) that a obeys the even parity-time reversal condition 

o-(r,0,0) = o-(-T,7r-M + 7r), (3-13) 

so that the standard covariant phase space symplectic structure is finite. Remark that this 
condition also completely removes the freedom of performing logarithmic translations since H 
are (parity-time) odd functions on the hyperboloid. In this work we will relax the conditions 
(3.11), (3.12) and (3.13) but impose that k ab is traceless and divergence-free. 

To derive the conserved charges associated with rotations and boosts, the specification of 
the fall-off conditions (3.1) are not sufficient. In this work, we will restrict ourselves to the phase 
space where the metric can be put in the following Beig-Schmidt form up to second order 

ds 2 = (l + j + ^+o(p~ 2 ^jdp 2 + o(p- 1 )dpdx a 

dx a dx h . (3.14) 




The logarithmic term i ab is necessary and sufficient in order that the equations of motion have 
a solution in general, up to second order in the expansion, when a and k ab do not obey parity 
conditions (see Section 4.3). The set of metrics (3.14) is a consistent truncation of the phase 
space (3.1), which is preserved by Poincare transformations, logarithmic translations and super- 
translations. Note that subleading terms in those transformations are constrained such that the 
form (3.14) is preserved. 
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3.2 Boundary conditions in Hamiltonian formalism 

It is generally accepted that asymptotically flat spacetimes in 3+1 formalism should be defined 
using the following fall-off conditions on the three-dimensional metric 3 gij and its conjugate 
momentum n" 13 

%j = 8 l3 +0(r- 1 ), d k 3 9ij = 0(r- 2 ), (3.15) 

n v = 0{r~ 2 ), 

on the initial data surface S defined as t = 0. Here Latin indices denote spatial components: 
x l = x 1 ,^ 2 ,^ 3 in Cartesian coordinates and r 2 = XiX i . We will use capital Latin indices 
A, B, . . . to denote both normal (_L) and spatial components and greek letters l, £, £ as indices 
on the two-sphere in spherical coordinates. 

Following [1], ten asymptotic values for the lapse and shift functions N A , M AB = M^ AB \ 
are introduced as additional canonical variables (together with their conjugate momenta) on the 
same footing as 3 gij, ix %3 . They specify the asymptotic location of the initial time Cauchy surface 
£ at spatial infinity through the asymptotic values of the lapse iV and shift functions N l 

N = M i ± x i + N 00 + O(r~ 1 ), 8iN = + 0(r- 2 ), 
N* = M/x j + N l OQ + Oir- 1 ), djN' 1 = Mf + 0( r - 2 ), (3.16) 

where 

M i:j = Mi k 5 kj = —Mji, M i ± = -M i± = M ±i . (3.17) 
In the original Hamiltonian treatment [1], the leading asymptotic parts of the canonical fields 

%j = 5ij + \*g$ + o(r- 1 ), (3.18) 

ir ij = ^ {2)ij + o(r~ 2 ), (3.19) 

are restricted to obey the parity conditions 

3 gff(-n) = 3 g£\n), vr^-n) = -7r( 2 )^'(n), (3.20) 

where n = r _1 (x 1 , x 2 , x 3 ). These conditions have been shown to be necessary in order that the 
Hamiltonian defined from the canonical variables 3 gij, k %3 , M AB , N A is finite when asymptotic 
rotations and boosts are considered. The full class of transformations preserving the fall-off and 
parity conditions is larger than the Poincare group. Indeed, angle-dependent shifts of the lapse 
and shift functions 

N 1 = S ± {n) + o{r°), N* = 5*(n) + o(r°), (3.21) 

are allowed transformations once they are restricted to be parity odd 

S ± (-n) = -S ± {n), S i {-n) = -S i (n). (3.22) 

We will refer to 5 ,_L (n) as temporal supertranslations, S r (n) as radial supertranslations and 
supertranslations tangent to r constant surfaces S^(n) as angular supertranslations. 
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In this work, we relax the parity conditions (3.20). Once parity conditions are relaxed, the 
fall-off conditions are preserved by both parity even and parity odd supertranslations (3.21), i.e. 
conditions (3.22) are no longer required. In addition, the fall-off conditions on the asymptotic 
fields 3 gij and are preserved under the so-called logarithmic translations which are associated 
to the lapse and shift functions 

N 1 - = log r K- 1 +o(r°), (3.23) 
N" 1 = log riT + o(r°), (3.24) 

where K 1 - and K % are constants. 

In order to derive the explicit formulae for the conserved Hamiltonian charges associated 
with rotations and boosts, it is required to specify the subleading components of the fields in 
(3.18)-(3.19). For this purpose, we restrict ourselves to the truncation of the canonical phase 
space where the fields obey 

% - «, + I^> + 4%» +*-*), (3.25) 
*« = + + L W«+„( r -3). (3.26) 

The logarithmic branch is necessary in order that the Hamiltonian and momentum constraints 
admit solutions in general when parity conditions (3.20) do not hold, as already noticed e.g. in 
Appendix B of [16]. 

A comparison of boundary conditions between 3+1 cylindrical and four-dimensional hyper- 
bolic representations can be made following [15, 16]. We detail this comparison in Appendix 
A. Note that Poincare and logarithmic translations can be mapped between 3+1 and covariant 
formalisms. There is also a one-to-one mapping between canonical temporal and radial super- 
translations and covariant supertranslations verifying (3.7). Covariant angular supertranslations 
generate mixed components g pa in hyperbolic asymptotic coordinates and are usually fixed in 
the Beig-Schmidt expansion [35]. As our approach is based on the existence of a variational prin- 
ciple, we only consider covariant boundary conditions transposed to the Hamiltonian formalism. 
We will thus relax conditions (3.20) and (3.22) but impose additional boundary conditions on 
the fields such that the angular supertranslations are fixed. These last conditions impose that 
the a priori generic form of the fields present in (3.25)-(3.26) is actually fixed in terms of their 
covariant counterparts as detailed in Appendix A. It would be interesting, however out of the 
scope of this paper, to include mixed terms g pa in the Beig-Schmidt expansion and allow for 
angular supertranslations along the lines presented in the next section 4 . 



4 Lagrangian dynamics in the hyperbolic representation 
4.1 Action principle 

In this section, we define the action principle for asymptotically flat spacetimes using hyperbolic 
temporal and spatial cutoffs. We introduce a finite hyperbolic cut-off p = A in order to regulate 

4 In the presence of mixed terms g pa , there might be a distinction between the bulk covariant phase space 
symplectic structure defined from the action and the one defined from the equations of motion, see definitions of 
lu and W in Section 4.4. One would then need to prescribe which one is the bulk symplectic structure, see [37] 
for an example where such a prescription plays an important role. 
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spatial infinity. We denote as % the unit hyperboloid defined at the cut-off p = A. We also limit 
the domain where the variational principle is defined between an initial and a final hyperbolic 
spacelike hypersurface that we denote as E±. Such temporal slices are relatively boosted with 
respect to each other close to spatial infinity. The spheres lying at the intersection of the 
boundary hyperboloid H with the hypersurfaces S± are denoted as S± and are defined at 
hyperboloid times r = t±, see Figure 1. 




Figure 1: The variational principle is defined in the spacetime delimited by initial and final 
hyperbolic temporal slices £± and the hyperbolic radial cut-off %. 

The notation we will use through Section 4 is as follows. We will denote the volume form as 

^- p x)^.., v = p , (n 1 _ p) , em -M»^ P+1 A • • • A dx^, (4.27) 

and we also have dg^ = g m g v ^ ] 5g a p. The covariant derivative associated to g pv is V^. Two- 
forms are defined as k = ^k pu dx^ A dx u = k^ u (d 2 x) pu implying that dk = ^Vak^dx' 7 A 
dx^ A dx v = V ' v k^ h '{d 3 x) p . For Stokes' theorem to hold, we have set conventions that imply 
^Tp6<f> = £p6<f> = — €t6»</> = £e<j>- This also tells us that (d 3 x) T = —dp(d 2 x) T , (d 3 x) p = —dr{d 2 x) p 
and (d 2 x) T = — (d 2 x) p = 2(d 2 x) pr = —d 2 S where d 2 S = ^e^dx 1 * A dx L . 

Using a regulation with hyperbolic cutoffs, the action for asymptotically flat spacetimes has 
the form 

S = ^Tnf d 4 x^R + Sx ± +S n , (4.28) 
1d7tG J m 

where Sy, ± are boundary terms at the future and past boundaries and Su are boundary terms 
at the radial boundary. It is interesting to first concentrate our attention on the past and future 
boundaries that are usually not treated in the formulation of a variational principle with zero 
cosmological constant. Upon varying the action, the boundary terms on "E± can be written as 

5S\x ± = ^ (d 3 x) p e»[5g] + 55 s± + B s± , (4.29) 

where 

e»[5g](d 3 x) p = ^g (g^6Y^ - g^5T%) (d 3 x) p (4.30) 
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is the presymplectic form obtained from 5{^—gR) = —\^—gG fiV Sg IJ , ll + M G^. The terms Bs ± 
defined on the boundary spheres S± are obtained from collecting the boundary terms in the 
variation of the action 5% with respect to boundary fields. Using the asymptotic expansion of 
the fields and f- s (d' i x) a Lo a = — J p ~ dp J s {d 2 x) a uj a , we obtain that the first term in equation 
(4.29) contains a linearly divergent term, a logarithmic divergent term and a finite term J-j:± 
that can be written as 

5S\v ± = J d 2 S^h^n a UaV a da + h bc V a 5k b )j + 6 (<S S± + ATZ S± ) 

+^ ± + f3 s± , (4.31) 

where TZs± ~ fg, d 2 Sy / —h^a. The finite term Fy, ± as well as the linearly divergent term 
ASlZs ± are expected to be exactly canceled by an appropriate choice of counterterm 5s ± after 
an appropriate choice of boundary conditions at future and past times has been made. The 
derivation of the boundary conditions and boundary terms at £± would require a careful analysis 
that we will not perform here. Now, the logarithmic divergent term at £± cannot be canceled 
by a local algebraic expression of the boundary fields alone. Indeed, there is no logarithmic 
branch in the asymptotic expansion of the metric (3.1) at leading enough order in p. We will 
discuss two mechanisms to set that divergence to zero. 

One mechanism to cancel the logarithmic divergence is to impose the following parity con- 
ditions 

<r(T,e,4>) = * ff ff(-T,7r-0,0 + 7r), (4.32) 
k a b(T,9,(/>) = s k k ab (-T,ir -6,(p + ir) . (4.33) 

and choose the spatial boundary counter-term such that Bs ± is finite or zero. Here, s a and 
Sfc are two signs which define the phase space with parity conditions. The Hamiltonian parity 
conditions imposed in [1] amount to s a = Sk = +1 5 . 

Let us now show that there is another mechanism for canceling the logarithmic divergence 
without assuming that a or k ab obey parity conditions. We simply propose to choose the spatial 
boundary counter-term to be 

S H = S% rt + (s^ + S«) , (4.34) 

where and S<® are the actions 

= I d 3 xv^hW (i(7(D + 3>) +\ I {d 2 x) a aV a a- l -j (d 2 x) a aV a a 

= [ d 3 xV-/i(°) ( - 1 V a aV a a + 3 a 2 ) + / (d 2 x) a aV a a - [ (d 2 x) a aV a a, (4.35) 
Jh v 2 2 / J s+ J s _ 

S {k) = f^x^/^)^k ab (U-?,)k ab ) + ^f s {d 2 x) a k cd V a k cd - 

= jf d?xsfW) (~ o6 5«) + (d 2 x) a k cd V a k cd - 1 jT (d 2 x) a k cd V a k cd (4.36) 

5 A necessary and sufficient condition to allow the Schwarzschild black hole is s CT = +1. The choice Sfc = +1 is 
necessary in order to disallow parity-odd supertranslations and therefore non-trivial supertranslation charges, as 
can be shown from the explicit expression (4.87) below. Together, these parity conditions are equivalent to the 
parity conditions of Hamiltonian fields imposed in [f ] , see dictionary in Appendix A. 
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for the scalar a and the traceless, divergence-free tensor k ab defined on the hyperboloid. Also, 
S^ rt is a counter-term action built algebraically from the boundary fields (and that therefore 
contain no logarithmic divergence in Bg ± ) and that is chosen to cancel the boundary terms at 
spatial infinity, e.g. the Mann-Marolf counterterm [3]. We will discuss that counter-term action 
below. 

The variations of the actions and S^ k ' reduce on-shell to a boundary term at S± because 
the equations of motion for SW and are precisely the equations obeyed by a and k ab obtained 
from Einstein's equations, see Section 4.2. Now, the sum of logarithmic divergences of the bulk 
term and boundary counter-term at S± precisely cancel for the choice of coefficients in (4.34) 
and the choice of boundary terms in the actions (4.35)-(4.36). 

The action (4.34) explicitly breaks translation, log translation and supertranslation invari- 
ance but does not break Lorentz invariance. The presence of logarithmic counter-terms is remi- 
niscent of the Weyl anomaly [38, 39] in the holographic renormalization of anti-de Sitter space- 
times in odd spacetime dimensions [21]. We will refer to the action 

A = — *— [ tfW-/i (0) fcr(D + 3)<r + 1 k ab {U - 3)k ab ) + (boundary terms) (4.37) 
8itG J% V 8 / 

as the (super/log-)translation anomaly. The anomaly is invariant under all symmetries that 
are broken. Indeed, translations do not act on the fields a and k ab . Logarithmic translations 
and supertranslations act as (3.4) and (3.5) but the anomaly is invariant up to boundary terms 
at timelike boundaries (that we neglect for this argument). Therefore, the Noether charges of 
the anomaly associated with the (super/log)-translation symmetries represent the algebra of 
(super/log)-translations. The Wess-Zumino consistency conditions [40] are therefore obeyed. 
Even though the anomaly is zero on-shell for all metrics obeying the boundary conditions, it 
affects the dynamics mainly because its symplectic structure is non-zero on-shell, as we will 
discuss in Section 4.4. Since no holographic model for asymptotically flat spacetimes is known, 
we cannot unfortunately try to match the flat spacetime anomaly to a QFT model. 

When parity conditions hold, no logarithmic divergent term appears in the variation of 
the Einstein-Hilbert action. There is therefore no requirement to add an anomaly term in the 
Einstein-Hilbert action. If one insists in defining an unique action principle whether or not parity 
conditions hold, the phase space of asymptotically flat gravity where parity conditions do hold 
also has a non-vanishing anomaly in the action. The variation of the anomaly does not lead to 
any logarithmic divergent boundary term because of parity conditions, which is consistent with 
the fact that there is no divergent term to cancel. The anomaly therefore does not contribute 
to the symplectic structure and the conserved charges, as we will detail below. More generally, 
the anomaly does not participate in any on-shell dynamics. For all purposes in classical gravity, 
one can therefore ignore the anomaly. The anomaly might however have a role in a quantum 
path integral involving off-shell configurations. Since this work is purely classical, we will simply 
discard the anomaly when the phase space is restricted with parity conditions. 

Let us now discuss the remaining spatial counterterm S^ rt . Following earlier work [6, 8, 9, 
10, 11], it was proposed in [3] to define 5^ as 

S^ = -±-f d'xV^h(K-k), (4.38) 

07TG J p=A 

where K is the trace of the extrinsic curvature, K = h ab K ab and K ab is defined implicitly by 

iz ab = k ab k - k a c k cb , (4.39) 
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where 7Z a b is the Ricci tensor of the boundary metric h a b- This equation being quadratic in 
K a b, it admits more than one solution for K a b- The prescription of [3] consists of choosing the 
solution that asymptotes to the extrinsic curvature of the boundary of Minkowski space as A is 
taken to infinity. It can then be shown that the variation of the action at the spatial hyperbolic 
boundary is equal on-shell to 

SS\ P=A = ^ j tfx^nm ab 5k ab , (4.40) 

where E^) = — cr & — cr/r? is the first order electric part of the Weyl tensor. The variation of the 
anomaly action (4.37) is zero on-shell (up to the crucial boundary term Bs ± treated earlier) and 
therefore does not contribute on the boundary hyperboloid. Here, Einstein's equations imply 
that h^ b is locally the metric on the hyperboloid and therefore we set SnfJ = by fixing the 
boundary metric to be the unit hyperboloid. 

Now, let us observe that the condition 5k a b = imposed in [3] can be relaxed. Using 
integrations by parts, the boundary conditions (3.6) imply that 

5S\ P=A = — l — [ d 3 xV-hWa5k = 0. (4.41) 

As already mentioned, we expect that one can constraint the initial and final states and add 
appropriate boundary terms at past and future times in order to cancel the boundary terms 
at £± in the variation of the action 5S. Therefore, up to that assumption, the action (4.28) 
is a valid variational principle for all k a b obeying (3.6) without necessarily imposing parity 
conditions. Note that the trace k is conjugated to the mass aspect ratio a in the variational 
principle, which is reminiscent of the AdS variational principle where the boundary metric is 
conjugated to the stress-energy tensor [24, 25] 6 . 

If one insists in imposing parity conditions on a and fc a &, the covariant phase space defined 
by our boundary conditions is still more general than the one considered in [2, 34] since here 
we do not impose k a b = but only require that T> a k a b = k% = 0, which implies that the 
first order part of the magnetic Weyl tensor B^} = ^e a cd, D c kdb can be fluctuating and non- 
vanishing 7 . In our phase space with parity conditions, the parity even supertranslations fulfilling 
(3.7) are allowed while logarithmic translations are forbidden. As we will see herebelow in 
Section 4.5, the conserved charges associated with those allowed supertranslations are vanishing 
while the Poincare generators are non- vanishing. Therefore, the phase space defined by our 
boundary conditions (3.14) supplemented by parity conditions is a consistent phase space where 
the asymptotic symmetry group is the Poincare group. The covariant phase space defined by 
(3.14) with parity conditions is a subset of the phase space defined in 3+1 formalism by Regge 
and Teitelboim in [1] where angular supertranslations have been fixed by requiring g pa = °(p~ 2 )- 

In the following, we will derive the dynamics of the phase space defined from our bound- 
ary conditions from the action principle. Since we regularized all divergences the symplectic 
structure and conserved charges will be finite when all boundary contributions are taken into 
account. 

6 We thank K. Skenderis for pointing this out. 

7 (1) 

Note that the four lowest harmonics in B ab are zero because we impose that k a t, is regular, see Lemma 1 in 
Appendix D. 
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4.2 Equations of motion 

Using the 3+1 split of Einstein's equations along the hyperbolic radial coordinate p and using 
our expression for the metric (3.14), we can expand the equations up to second order in the p 
expansion. This analysis is performed in Appendix C. By convention, we refer to the equations 
involving only h *~) as the zeroth order equations, the ones linear in k ab , a as the first order 

(2) 

equations, and the ones linear in i ab and h , as second order equations. The zeroth and first 
order equations were already obtained in full generality in [35]. As already said before, the 
zeroth order equations imply that h$ is locally the metric of the unit hyperboloid. At first 
order, setting k ab to be traceless and divergenceless, we obtain 

(□ + 3)<r = 0, (D-3)k ab = 0. (4.42) 

One can rewrite elegantly those first order equations in terms of the first order electric and 
magnetic parts of the Weyl tensor 

= -V a V b a-h<$a, B^ = \e a cd V c k db . (4.43) 
By construction, these two tensors enjoy the following properties 

h ( )a bB (l) = Q > E W = > V[cE W = . (444) 

1) 

■lb] 



The boundary conditions (3.6) imply that B^ = 0. The equations of motion at first order 
(4.42) are then equivalent to 

h ( ° )ab E£! = , V [C BU = 0. (4.45) 

At second order, the equations derived in the literature do not include k ab and i ab . A 
complete derivation is provided in Appendix C and the resulting equations are given from 
(C.167) to (C.172). Instead of giving their explicit cumbersome expressions, and following the 
works presented in [36] and [12], these equations can be reformulated in an interesting and 
compact way by introducing a pair of conjugate symmetric, divergence-free and traceless (SDT) 
tensors V ab and W ab . They are defined as 

Vab = -h^ + hab + Q^b, (4.46) 



U,„ ee e a cd D c (h^-\i db + Qf b ) (LIT) 



where QY, and Q^ b are non-linear terms defined in appendix in (C.187)-(C.182). The tensors 
are conjugate in the sense that they obey the following duality properties 

W a6 = -(curiy) o6 + (curlre) (a6) , V ab = (cuAW) ab - 2i ab , (4.48) 

where 

- = 4r' 1 +14r , +>•^+4• M, (<>■*» 
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is a sum of symmetric and divergence-free (SD) tensors. All SD tensors built out of quadratic 
terms in a, k ab and their derivatives can be completely classified. This classification is performed 
in Appendix B. The particular tensors appearing in K ab (4.49) and some of their important 
properties are described at length in Appendix B. 

Using the SDT tensor W ab , it is shown in Appendix C.3 that the equations of motion can 
be written as 

W£ = v b w ab = o, 

(n-2)W ab = cml{2i + K) {ab) , (4.50) 
i a a = V b i ab = 0, 
(□-2)t o6 = 0. (4.51) 

Using the curl operator and the definition j ab = —(curli) ab , one can derive an equivalent form 
of those equations in terms of V ab as 

v a a = v b v ab = o, 

(D-2)V ab = curl {-2j + curl (k)U, (4.52) 

f a = V h j ab = ^ 
(D-2)j ab = 0. (4.53) 

Since the curl of the latter set of equations lead to (4.50)-(4.51), the two sets are equivalent. 

(2) 

Once the set of equations (4.50)-(4.51) is solved, one can reconstruct h J from the definitions 
(4.46) or (4.47). Einstein's equations at second order are therefore equivalent to either set of 
the above systems of equations. 

4.3 Linearization stability constraints 

For any SDT tensor T ab and any rotation or boost vector £? Q s, such that ^( a C(o)6) = an d 
(□ + 2)^ 0) = 0, we get 

(□ - 2)T ab ^ 0) = 2 V a ($ Q) V [a T b]c + T c[a V b] e {0) ) ■ (4.54) 
Therefore, the following integral 

! d 2 S (□ - 2)T ab efo n b = (4.55) 
J s 

vanishes as the r.h.s. of (4.54) reduces to a total divergence on the two-sphere. As pointed out 
in [36], necessary conditions for the equations of motion (4.50) to admit solutions are that the 
r.h.s. of (4.50) contracted with and integrated on the sphere vanishes 

J d 2 S (□ - 2)W ab £f 0) n b = J d 2 S curl (2i + k) (o6) ^ 0) n b = . (4.56) 

Such conditions were referred to as integrability conditions in [36]. They can be recognized as 
linearization stability constraints [41, 42]. 
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Using the fact that charges constructed with the curl of an SD tensor (such as n a b or i a &) 
can be mapped to charges associated with the SD tensor itself, as proved in the appendix B of 
[12], these conditions are equivalent to 

/ d 2 Si ab i a m n h = -\f d 2 S Kab £f 0) n b . (4.57) 
J S J s 

Now, using the properties of K a b detailed in Appendix B, it turns out that there exists a unique 
equivalence class of tensor structures quadratic in a, respectively in k a b, which might evaluate 
to a non-zero result in the r.h.s of (4.57). For the mixed (cr,k) SD tensors, we show that none 
of them could contribute. 

Remarkably, the two equivalence classes of charges quadratic in either a or k a \, can be 
recognized as the Noether charges associated with the two free actions (4.35)-(4.36). This 
provides an independent argument that the actions (4.35) and (4.36) play an important role in 
the dynamics of asymptotically flat spacetimes. More precisely, using the Noether charges of 
the actions (4.35)-(4.36) as representatives of the two classes of conserved charges that one can 
build out of quadratic terms in a and k^, we can rewrite the linearization stability constraints 
as 

d 2 Si ah e (0) n b = 2 [d 2 S (t£> + r«W 0) n 6 , (4.58) 



where 



_ 2 5L^ _ 1 [ a ,a,I] 1 [a,ff,II\ 



{k) L 6LW _ 1 [k,k,II] l y (2) 



are the stress-tensors of the actions and S^ k ' given in (4.35)-(4.36). The SD tensors ap- 
pearing on the r.h.s of (4.59) and (4.60) are defined in Appendix B. 



4.4 Symplectic structure 

In this section we first review the standard covariant phase space symplectic structure [43, 44, 45]. 
We then fix the well-known ambiguity [18, 19] in its definition using the boundary terms in the 
action principle (4.34). Effectively, we will fix the boundary terms in the symplectic structure 
in such a way that the logarithmic divergences, present when parity conditions are not imposed, 
cancel. Our procedure then amounts to the prescription already proposed in [20]. The final 
symplectic structure that we define has the form 

tt[5 1 g,5 2 g] = fWk[<5iS, S 2 g] + Uct.^io", hk; S 2 ct, S 2 k], (4.61) 

where Qbulk is the standard bulk symplectic structure and £l c .t. is precisely the boundary sym- 
plectic structure that can be derived from the boundary action in (4.34). 

The symplectic structure is a phase space 2-form defined as an integral over a Cauchy slice 
of a spacetime 3-form that we refer to as the integrand for the symplectic structure. The 
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Ashtekar-Crnkovic-Witten-Lee-Wald integrand [43, 44, 45] for the covariant phase space sym- 
plectic structure is given by 

w[8ig, hg\ = ^^{d z x)^g( K 5 1 g a ^5 2 g a p + 5 ig V a 5 2 g^ + 6 ig ^7 a S 2 g 

-5 ig V^ 2 g - 26 l9a pV a 5 2 g^ - (1 o 2)) , (4.62) 

where 5\g^ u , 5 2 g fiu are perturbations around a general asymptotically flat spacetime g^. This 
integrand is obtained by varying a second time the boundary term &[8g] obtained after a varia- 
tion of the Einstein-Hilbert Lagrangian as Lo[S±g, S 2 g] = 5i&[5 2 g] — 5 2 ®[5ig\. The integrand for 
the symplectic structure, obtained in [46] by acting with a contracting homotopy on Einstein's 
equations contracted with 5g^ u , is given instead by 

W[Sig, 5 2 g] = v^S P» vaM& V »6 2 g j5 - 5 2 g a/3 V u 5 2 g l5 

= ^(d 3 x)^(Si9 a ^5 2 g a p + 5 ig V a 5 2 gZ + 6 ig ^V a 5 2 g 

-hgV^zg - 5 iga pV a 5 2 g^ - 6 ig ^5 2 g a p - (1 O 2)) , (4.63) 

where 

_ gfj,v g-y(a @)6 _|_ fi(-y g 5)u g aP _|_ fi(a 0)v g -y8 

_ gl *v g af} gl 5 _ g K-y g S)(a g P> _ g K<* g P)(n g &> , (4. 6 4) 

The integrands of the two symplectic structures defined above differ by the boundary term 

W[6ig,S 2 g]-cj[6ig,S 2 g} = -J-^a^V^V v {d ig ^6 2 g^ - (/x -> i/)), (4.65) 

which vanishes on constant p and r surfaces when the metric and perturbations are expanded 
in the Beig-Schmidt expansion with 

g P a = 0, 5g pa = . (4.66) 

We can therefore use interchangeably oj and W in what follows. The integrand for the symplectic 
structure derived from Einstein's equations contracted with the Lie derivative of the metric is a 
boundary term 

W{6g,£^g) = dk^[Sg;g], (4.67) 

where the last equality has to be understood up to Einstein's equations of motion for the metric 
and the linearized perturbations. The boundary term k^[Sg;g] is exactly given by the Abbott- 
Deser expression [31] for the surface charge for a linear perturbation Sg pu around a solution 
g^v 

' ;(d 2 x) flu ^U u (^Sg - V a 6f*) + ^Vi/" 
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+ l -8gV v e + + W'V^ - ( M v) ) . (4.68) 
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Let us now discuss conservation and finiteness of the symplectic structure. On the one 
hand, using the expansion (3.1) at first order, the symplectic structure integrand evaluated on 
a hyper surf ace p = constant gives 

W[8 ig ,6 2 g]\ fixedp = o(p ), (4.69) 

which shows that the symplectic structure is conserved. On the other hand, the integrand for 
the symplectic structure evaluated on a Cauchy slice £ asymptotic to a constant r hypersurface 
reads 

W[6 19 , 5 2 g]\x = £^ 3 ^ V C /^ (fioVfo ~ ~/ ce 8ik c d 5 2 B« - (1 2)) + o^ 1 ), 

(4.70) 

which is consistent with the earlier expression (4.31). The bulk symplectic structure 

n Mk [S ig , S 2 g] = £ W[S 19 , 5 2 g] (4.71) 

is therefore logarithmic divergent for generic a and k a f,, in accordance with the variational 
principle being ill-defined on future and past boundaries £± when parity conditions are not 
imposed. 

We have seen in Section 4.1 that there is a mechanism to cancel the logarithmic divergence 
in the action without imposing parity conditions. The logarithmic divergences are canceled once 
the Einstein-Hilbert action is supplemented by a logarithmic counter-term, as in (4.34). Now, 
even though the logarithmic counter-term action is zero on-shell (up to boundary terms), it has 
a non-vanishing boundary contribution to the symplectic structure. Indeed, the free actions for 
a and fc a j, introduced in (4.35)-(4.36) are zero on-shell but their symplectic structures (defined 
with the same conventions as in the bulk) are the Klein-Gordon norm and the symplectic norm 
between two traceless transverse fields given by the integral of 

0J {a) [5 1( j,8 2 a] = (d 2 x) a V-h(°) (8 1 a& 1 5 2 a - (1 O 2)) , (4.72) 

w (fc) [<5iM 2 fc] = ^^^(^^^^-(l^^, (4.73) 

over the sphere which is generally non- vanishing 8 . The total symplectic structure is then defined 
as announced in (4.61) with 

n c . t . [Skt, S t k, 8 2 a, 8 2 k] = ^ ! [6 lt r, 5 2 a] + [Sik, 8 2 k]). (4.74) 

J s 

The resulting prescription for fixing the boundary terms in the symplectic structure left unfixed 
in [45, 18, 19] amounts to the prescription argued in [20] to fix the boundary terms in the 
symplectic structure using the symplectic structure of the boundary terms of the action. We see 
that we can now justify that prescription by the existence of a variational principle when past 
and future boundaries are taken into account. 

8 The counter-term / _ A d s x\/—h(K — K) has not the form of an off-shell action for the boundary fields and 
therefore it does not define a boundary symplectic structure. 
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4.5 Covariant phase space charges 

In this section, we compute the covariant phase space charges, paying particular attention to 
the contributions coming from the boundary symplectic structure (4.74). Since the symplectic 
structure is finite, the conserved charges are also finite. Along the way, we review how the 
infinitesimal charge coming from the bulk symplectic structure (which is not finite in itself 
when parity conditions are not imposed) can be expressed equivalently in the form derived by 
Abbott-Deser [31], Iyer- Wald [18, 19] or Barnich-Brandt [33]. After some general considerations, 
we respectively consider the conserved charges associated with translations, supertranslations, 
Lorentz charges and logarithmic translations. When parity conditions are imposed, in addition 
to our boundary conditions, the resulting expressions for the charges reduce to the ones derived 
in previous works on asymptotically flat spacetimes in Lagrangian framework [29, 5, 30, 31, 32]. 
Indeed, as we have already stated, the only difference with those works, when parity conditions 
are imposed, is that we have allowed for temporal and radial supertranslations whose associated 
charges are zero. 

The covariant phase space infinitesimal charges associated with a diffeomorphism tangent to 
the phase space are defined from the symplectic structure as 

fQt\g] = n[6g,£tg]. (4.75) 

Here, we use the symbol fQ^[g] to remind the reader that the infinitesimal charge between the 
solution g and g + 5g can be considered as a one-form in field space which is not necessarily 
exact. When $Qz\g\ is an exact form in phase space, as it will turn out to be the case for each 
diffeomorphism that we consider, the charges can be defined. We denote the integrated charge 
as Q^[g; g] (we have 5Q^[g; g] = $Qt\g\) and fix the integration constant so that Q(.[g; g] = for 
Minkowski spacetime g. 

Using the definition (4.61) and the property (4.67) of the bulk symplectic structure, the 
charge one- form SQ^[g] can be written as a surface integral 

$QM = J h[5g;g] + ^ J (u,&[6*, Sp] + U W[6k, 8^]) , (4.76) 

evaluated on the sphere S at constant time t and at p = A. Here, 5^a, S^k a b are variations 
of the first order fields induced by the Lie derivative of the metric along £. The Abbott-Deser 
expression can be rewritten in the alternative form 

kdSg;g] = -«f + ^(<^Vv^(nv^- v^") + vac) -£[£ e5 ,^],(4.77) 

where 

kf[ 9 ] = ^(^^(VT - V"£") , (4.78) 

is the Komar term and 

E[Ct9, 8g] = x ^(d 2 x)^^-g( i -g^5g a pg^C i g 1<T g^ - ( M O vj) , (4.79) 

is a term linear in the Killing equations that might not vanish in general for asymptotic sym- 
metries. Here 5 acts on the metric and on the asymptotic Killing vector, as might depend 
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on the metric. Remark that in (4.77), the first term is the exact variation of the Komar term 
and the third term is zero when evaluated on constant p and r surfaces since 



9 P a = 0, 5g pa = 0, C^g pa = 0. (4.80) 

Therefore the bulk surface charge one- form (4.77) is given by the generalization of the Iyer-Wald 
expression when the asymptotic Killing vector is allowed to depend on the metric 

kt[6g;g] = -5kf[g] + ^ d \)^^g(e {^8g - V a 5g^) + V^f)- (4.81) 

The explicit evaluation of the total charge one-form (4.76) for each asymptotic Killing vector 
(translations, supertranslations, rotations, boosts and logarithmic translations) is a straightfor- 
ward exercise that we do not reproduce here in details. In the following, we simply present the 
key steps in the evaluation of the charges and state the final results. 

Translations 

For translations, it is relatively simple to see that 

[Sg; g] = % [5g; g] , <%cr = <% k ab = 0. (4.82) 

As a result, the boundary contributions to the charges vanish and the charges can be defined 
using linearized perturbations of asymptotic fields around Minkowski spacetime 

Qd9i9\= / fy[h;g], h pu = g pu - g^ u . (4.83) 
Js 

Effectively, one can therefore derive the charges associated with translations using the bulk 
linearized theory and obtain the well-known expression for the four momenta 

QV)\9i5] = -g^ / d 2 SE^n a V b ( M , (4.84) 
J S 

where the four scalars Q^), \x = 0, 1, 2, 3 are the four solutions of T> a T>f,(/ fJ \ + h^Q^ = 0. The 
vector d/dt in the 3+1 asymptotic spherical coordinate system (t,r,9,(f>) corresponds to the 
translation Q ) = — sinhr. One can check [14] that the charge Q(o) f° r the Kerr black hole of 
mass m is +m after identifying a = Gm cosh 2r sech r, as it should. 

Supertranslations 

For supertranslations, we start by showing that 

k([6g;g] = k^[5g;g], <%cr = 0, 5^k ab = 2(u ab + h^u). (4.85) 

The bulk charges can again be obtained using the linearized theory. Now, it is important to 
note that the boundary counter-terms do not contribute to the supertranslation charges 

[6a, Sp] =0, f u, k) [5k, 5 ( k]=0. (4.86) 
s Js 

The first equality is trivial. The second one can be proven after using 6^B^ = 0, B^) = 
—a ab — h^a, where (□ + 3)5" = as a consequence of Lemma 1 of [12] (see also [5]) and 
eventually the fact that n a e ace u^a e d = n a T> c I e ace (\oj d a e d — \a d ui e d — w5 e ) J is a boundary term. 
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The supertranslation charges are therefore precisely the ones which can be obtained from 
the bulk linearized theory. They read as 



Q(uj) [gi 9} = ^ d 2 Sy^hWn a (a a u - au: a ) . (4.87) 
J s 

As we restrict ourselves to supertranslations satisfying (3.7), the charges are conserved since 
V a (a a uj — auj a ) = Haul — aduj = 0. Note that if we choose ui to be a translation, the charges 
indeed reduce to the earlier expression (4.84) upon using 2u:ai ) — 2uji } a = —E b a u a — 2'D a (uj[ a a b j ). 
When parity even conditions are imposed on both a and k ab , uj also has to be parity even 
and supertranslation charges identically vanish. This does not apply to translations u = Q^), 
which are of odd parity, as these are still allowed and are associated with non-trivial charges. 
The supertranslations are therefore associated with non- vanishing charges only when the parity 
conditions (4.32)-(4.33) with s a = Sk = 1 are not imposed as a boundary condition. 

Lorentz charges 

The derivation of the Lorentz charges is much more involved. First, the property of asymp- 
totic linearity 

k ( [6g;g] = k € [Sg;g], (4.88) 

does not hold for Lorentz transformations £ = — £(o) 9 • Moreover, the boundary fields vary as 

<%cr = -£ 5(0) cr, S^kab = -£^ 0) k ab , (4.89) 

and the boundary counter-term charge does contribute. The computation goes as follows. The 
last term in (4.81) is explicitly zero for a rotation since £/n) = £(o),rot is then tangent to the sphere 
and does not depend on the metric. This is the familiar result that the bulk charge associated 
with rotations is given by a Komar integral. However, in this case, the Komar integral contains 
a logarithmic divergence as well as a finite piece. For boosts £( ) = £(0), boosts the last term in 
(4.81) contributes and, therefore, the boost charge one-forms are not manifestly exact. It turns 
out to be manifestly exact only after we have substituted the Beig-Schmidt expansion in this 
expression. In intermediate expressions, a linearly divergent term in p appears but is canceled 
between the Komar integral and the second term in (4.81). The expression for the bulk charge 
eventually admits a logarithmic divergence and a finite piece as in the case of rotations. 
The logarithmically divergent piece of the Lorentz charges takes the form 

= / d 2 Sv^hWi ab ^ b + O(A ). (4.90) 

J s 

It can be expressed using the linearization stability constraints (4.58) as 

k ^ 9]9] = / ^^(^ + T %%)" b + 0(A°), (4.91) 

J s 

where and are the stress-tensors of the actions (4.35)-(4.36). Following our prescription, 
this divergence should be exactly canceled by the counter-term contributions to the total charge 

9 We follow here the usual sign conventions of [19]. 
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(4.76). We have checked that the divergence indeed cancels after using the following relationship 
between symplectic structures and Noether charges, 

u ia) (6*,£-t m a) = 5( K V^h^T^ ab C {0)b (d 2 x) a ) - d 2 sV^n a V b (2£ {0) [a a b] 5a)(4.92) 

u (k) (5k,£-t m k) = 5 (^)T« a6 e ( o) ^ 2 x)a) + d 2 S^) n ^ 6 P H , (4.93) 

where P a b = P\ a b] is an anti-symmetric tensor. 

Finally, the Lorentz charges are finite and have the following form 



*} S 

+hf b \8a 2 + a c a c - h ab k ab + k cd a cd )) ^ 0) n b . (4.94) 

The charges are also conserved as a consequence of the momentum equation (C.171). When 
parity conditions hold, the integral of all quadratic pieces vanish and the integral of i ab also 
vanishes as a consequence of the linearization stability constraints, see e.g., [14]. It then implies 
that asymptotic linearity (4.88) holds and the counter-terms trivially integrate to zero. The 
charges thus agree with the Abbott-Deser formula 



Q-^ [9; g]= fc-e {0) [9 -g\g]< (4.95) 

J S 

When a and k ab are not parity even, the charges (4.94) contain quadratic terms in the fields 
and cannot be obtained from the linearized theory alone. 

We showed in our previous work [12] that the Lorentz charges can always be written in two 
equivalent ways using the conserved tensors V ab and W ab in the restricted phase space where 
k a b = iab = 0. The tensor V ab can be related to the second order part of the boundary stress- 
tensor [14]. Let us remark here that, after using properties of integrals of tensor fields on the 
hyperboloid proven in Appendix B of [12] and the equations of motion, we can rewrite the 
rotation and boost charges (4.94) in the equivalent forms 

J ^ ~ 8^gJ s d2S V^(V ab + 2i ab )g om n b = ~^qJ s d 2 Sv^mWab^ o S t {i) n b , 
K ® ~ ^j/ S "f^ ) ^ b + 2i ^°°^ nh = 



where V a b, W a b are defined in (4.46)-(4.47). One can check [14] that for the Kerr black hole, one 
gets the standard result = +ma for £(o),rot = J^- 
Logarithmic translations 

Logarithmic translations are allowed asymptotic transformations. They modify the first 
order fields as 



5^ = H, 5t.k ab = 0. (4.97) 

Since logarithmic translations do transform a, the boundary terms in the symplectic structure 
might play a role. Using integrations by parts one finds 

Ha b - H b a = -\e^ H a - V a (H [a a b] ) (4.98) 
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and therefore 



^ Ju^SaM*) = ^ J/S^mE%H«n» (4.99) 

where we have discarded the total divergence term on the sphere. The bulk covariant phase 
space charge associated with logarithmic translations is given by 

! ^[5g;g] = - l ^5 f d 2 ; sV^ ~e£ ] <H a n b + ^--5 f d 2 S V^hW k ab n a H b . (4.100) 
J S J S J s 

We find that the two divergent contributions are opposite of each other and exactly cancel. The 
remaining finite part is trivially integrable. Logarithmic translations are therefore associated 
with the non-trivial charges 



Q(H) = I d 2 sV^hWk ab n a H b (4.101) 



IS 

which are conserved thanks to the property T> a (k ab H b ) = 0. 

In the restricted phase space where k ab = 0, logarithmic translations are associated with 
zero charges or equivalently they are degenerate directions of the symplectic structure. When 
parity conditions are imposed, logarithmic translations are not allowed transformations and the 
associated charges do not exist. The presence of non- vanishing conserved charges associated with 
logarithmic translations is therefore a particularity of the phase space without parity conditions 
and with k ab ^ 0. 



4.6 Algebra of conserved charges 

In the last section, we obtained explicit expressions for conserved charges associated with trans- 
lations, Lorentz transformations, logarithmic translations and supertranslations. We obtained 
that all asymptotic charges are non-trivial in general in our phase space. The set of infinitesi- 
mal diffeomorphisms form a Lie algebra defined from the commutator of generators. A natural 
question to ask is whether or not the algebra of translations, logarithmic translations, Lorentz 
transformations and supertranslations is represented with the associated conserved charges. 

General representation theorems are available [47, 46] but one quickly realizes that they do 
not take into account boundary contributions to the symplectic structure. These contributions 
can be dealt with as follows. Every diffeomorphism in the bulk spacetime induces a specific trans- 
formation of the boundary fields through the Beig-Schmidt asymptotic expansion that identifies 
boundary fields from bulk fields. Therefore, the Lie algebra of infinitesimal diffeomorphisms 
defined from the commutator of generators also induces a Lie algebra of transformations of the 
boundary fields. The Poisson bracket between two charges is then defined as 

{Qd9;9],Q?{9;g]} = -5tQe\g;g], (4.102) 

where the variation 5^ acts on the bulk fields as a Lie derivative and on the boundary fields as the 
transformation induced on the boundary fields from the Lie derivative of the bulk fields. It would 
be interesting to develop a general representation theorem which takes boundary contributions 
into account along the lines of [48]. In this work, however, we simply evaluate the Poisson 
bracket using the explicit expressions for the charges derived in the previous section and taking 
into account the boundary field transformations. 
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Under an asymptotic translation £ = uj(x)d p + o{p°) where T> a T> b uj + h ab uj = 0, the boundary 
fields transform as 

S^a = 0, <Uto = 0, (4.103) 
6„iab = 0, S UJ V ab = V c (E^u c ) + 2e cd{a B^ )c co d (4.104) 

where and are the first order electric and magnetic parts of the Weyl tensor while 
Lorentz transformations £ = — £(o) act on the boundary fields as a Lie derivative 

6 -^ ) a = £-t(. ) a > 5^ (0) k ab = £-^ 0) k ab , (4.105) 

&-£((,) *a6 = ^-5 {0 ) *'a6, ^(0)^6 = ^-^0)^6- (4.106) 
Logarithmic translations act as 

fea = if, 5U = 0, 5 H i ab = -V c (E^H c )-2e mn {a B { b ] ) m H n (4.107) 

foKft = -lv c (k ab H c ) + 2V c (E^H c ) + 8e m lB$ m H n (4.108) 
and supertranslations act as 

<5 w a = 0, <5 u fc«6 = 2wo6 + 8j ab = (4.109) 

^ife = fc c (a w 6) + ^a&W + W C (P c /c a 6 ~ ^(ah)c) 

+ (°" c ^ c (a6) - o^ab - 2aujh^ + W (a CT 6 ) + CT c(a ^ } + ((T O w)) . (4.110) 

After an explicit evaluation, we find that all asymptotic transformations are well-represented: 
the Poisson bracket is anti-symmetric and is isomorphic to the semi-direct product of the Lorentz 
algebra with the (super)-translation and logarithmic translation algebra. In particular, the 
Poisson bracket between Lorentz charges and (super)-translation charges is given by 

{Q- kov 2(u,)} = -{Qh, G-* (0) } = 2 (u /), (J = £-e (0) w . (4.111) 
and the Poisson bracket between Lorentz charges and log-translation charges is 

{2-e (0)! Q iH )} = -{Q(h), 2- 5{0) } = G(H'), ff' = ( 4 - 112 ) 

Logarithmic translations and supertranslations obey the algebra 

{2(o,), Q(H)} = -{Q(H), S(")} = 4^ / #SyfjWn*(H*u-Hu a ) (4.113) 

where the right-hand side depends on the generators but does not depend on the fields. In 
the harmonic decomposition of u> on the sphere, the Poisson bracket is zero for all harmonics 
I > 1 and is a Kronecker delta for the four lowest harmonics I < 1. The algebra of asymptotic 
conserved charges is isomorphic to the algebra of asymptotic symmetries. No non-trivial central 
extension of the algebra is present. 

In order to derive these results, we used Appendix B and the Appendix B of [12] to simplify 
intermediate expressions and we discarded boundary terms. We have also used the property 
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described in [5, 35] (see also Appendix D) that regularity of k a b implies that the four conserved 
NUT charges [49] 

T>W = ^g [ #SB$n*V>>t w , V = 0, 1, 2, 3, (4.114) 
J s 

are zero. The presence of NUT charges would require modifications of the Poisson bracket, 
see e.g. [50] where a symmetric topological extension of the algebra of bosonic supercharges is 
considered. The fact that we find an anti-symmetric Poisson bracket gives us confidence in our 
intermediate expressions. 

We obtained that all transformations: translations, logarithmic translations, Lorentz trans- 
formations and supertranslations are well-represented despite the (log/super)-translation anomaly. 
The fact that the Lorentz group is well-represented is not surprizing given that the cut-off needed 
to regularize the action, see (4.34), is invariant under asymptotic Lorentz transformations. Now, 
it is also important to take into account the shifts in the action when one changes the cut-off 
used to regulate the action. These shifts can be analyzed as follows. 

Under a change of cut-off A, the action will be shifted by a finite piece proportional to 
the anomaly action + given in (4.35)-(4.36). The conserved charges associated with 
the asymptotic Killing vector £ will then be shifted by the boundary Noether charges of the 
action + associated with the symmetry 5^. Using standard manipulations 5^L = dK^, 
5L = j^$4> + d@[6(f)], the boundary Noether charges are defined as = — O One 
then quickly sees that translations and supertranslations are associated with vanishing Noether 
charges J s d 2 SJ^ = while Noether charges associated with logarithmic translations are propor- 
tional to the four-momentum Q(fj\ and Noether charges associated with Lorentz transformations 

are given by the integral of J $(0) = 2(T^ ab + T^ ab )^ 0) (d 2 x) a where T^ ab and T^ ab are the 
stress-tensors of the actions (4.35)-(4.36). 

Therefore, under a change of regulator, the translations and supertranslation charges are 
invariant. Logarithmic translation charges get shifted with the four-momenta and the Lorentz 
charges get shifted as 

AG-*<P)M - I d 2 Sv^(T^+T%% )n b . (4.115) 
J s 

These shifts can be obtained similarly by varying the regulator directly into the expression 
for the charges (4.91)-(4.100) before the subtraction of divergences between the bulk and the 
boundary. Four-momenta and supertranslations are finite without needing a regulator. They 
are therefore manifestly unchanged by the regulator. 

5 Hamiltonian dynamics 

We have seen that parity conditions on the hyperboloid are not required in order to define a 
consistent phase space in the hyperbolic representation of spatial infinity. Moreover, we have seen 
that when these conditions are relaxed, charges associated with Lorentz rotations and boosts are 
non-linear functionals of the first order fields and logarithmic translations and supertranslations 
are associated with non-trivial charges. Both these characteristics are not shared with the 
standard treatment of Hamiltonian charges at spatial infinity [28, 1, 16]. There, parity conditions 
on the sphere are imposed in order that the rotation and Lorentz boost charges are finite. Also, 
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charges are linear functionals of the boundary fields, logarithmic translations are not allowed 
transformations and supertranslations are associated with trivial Hamiltonian generators. The 
purpose of this section is to resolve this tension by proposing how the results of [1, 16] can be 
accommodated to enlarge the phase space to fields which do not obey parity conditions. 

The canonical two-form on the canonical phase space used in the treatments of [28, 1, 16] is 
the bulk canonical two-form 

fi(<Ji 3 g, 5nr, 5 2 3 g, 5 2 tt) = £ d 3 x( y 5 1 7r mn 5 2 3 g mn - 5 2 T: mn 5^g mn ) (5.116) 

defined from the bulk canonical fields 3 g mn and 7r mn at the initial time surface E at t = 0. In 
the case of asymptotically flat spacetimes without parity conditions the bulk canonical two-form 
suffers from a logarithmic radial divergence. Using the boundary conditions (3.18)-(3.19) one 
can express the canonical two-form as 

n(S 1 3 g,5 1 n,6 2 s g,6 2 ir) = (finite) + ^ f d 2 S (5^ mn 5 2 3 g^ n - <5 2 tt« mn 5 1 3 5 « ) (5.117) 

J s 

where A is a large radial cut-off and S is the sphere at r = A. Now, in complete parallel to 
the Lagrangian treatment, we propose to modify the dynamics by adding a boundary term to 
the canonical form. We proceed by first writing the boundary actions (4.35)- (4.36) at t = in 
the 2+1 decomposition (the boundary metric becomes the real time line times the unit sphere). 
We then switch to the Hamiltonian formulation of the boundary action and propose to supple- 
ment the bulk canonical fields with the canonical fields of the boundary Hamiltonian. We then 
introduce counter-terms to the canonical form in order to minimally cancel the divergences in 
(5.117), in complete parallel to the Lagrangian prescription (4.74). The regulation breaks trans- 
lation, supertranslation and logarithmic translation invariance. We interpret this breakdown as 
a consequence of the translation anomaly in the action, which is manifest only when fields have 
both parities. 

Let us now discuss briefly the form of the Hamiltonian generators associated with asymptotic 
Poincare transformations and supertranslations. The Hamiltonian generators contain two parts: 
the part coming from the bulk canonical form and the counter-term contribution that cancels 
the logarithmic divergences. The surface charge derived from the bulk canonical form associated 
with a gauge parameter e A = (e -1 , e m ) is given by [1] 

[5 3 g, Sit; 3 g, vr] = G mno ^(e ± 5 3 g opln - ef n 5 3 g op ) + 2e°5( 3 g on ir mn ) - e m 5 3 g no 7T no (5.118) 

where 

G mno P = 1 ^ ( Y l ° Y p + 3 9 mp Y° - 2 3 g mn 3 g op ) (5.1 19) 

is the inverse De Witt supermetric. Now, one can readily obtain that this expression admits at 
most a logarithmic divergence when one uses our boundary conditions. The explicit asymptotic 
expansions provided in Appendix A are useful in deriving the explicit forms of the charges. The 
logarithmic divergence is then exactly canceled by the boundary counter-term. The resulting 
final expressions for the charges in Hamiltonian formalism can then be obtained by a straightfor- 
ward explicit evaluation. We will not provide them here. We only note that the four-momenta 
are given by the usual ADM formulae, while the charges associated with rotations and boosts 
contain non-linear contributions in the canonical fields. 
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6 Summary and Discussion 



We showed that the four-dimensional Einstein-Hilbert action admits a logarithmic divergence at 
spatial infinity when one regulates the action using hyperbolic temporal slices or, in other words, 
when one considers the action principle between initial and final time slices that are relatively 
boosted with respect to each other at spatial infinity. We described two alternative mechanisms 
to cancel the divergence: either one imposes parity conditions on the first order boundary fields 
or one adds to the action an anomalous counterterm that breaks asymptotic diffeomorphism 
invariance including asymptotic translation invariance. We will summarize what each choice of 
boundary conditions implies. 

If one insists in requiring asymptotic Poincare invariance, one is led to impose parity condi- 
tions. In that case, we built a bigger covariant phase space than the one previously considered 
[5, 2] while preserving all essential properties of asymptotic flatness such as the Poincare group 
as asymptotic symmetry group and the inexistence of other conserved charges at spatial infinity 
than the Poincare charges. The enhancement of the covariant phase space that we have con- 
structed with respect to the one of [5, 2] originates from the admissibility of fluctuations of the 
first order part of the Weyl tensor (while still imposing regularity of all fields and therefore no 
NUT-momenta) and from the presence of temporal and radial supertranslations that are pure 
gauge. We will further comment on this generalized phase space below. 

If one insists in building the largest possible phase space, one is led to relax the parity condi- 
tions. This set of asymptotically flat spacetimes has new qualitative features that we unraveled. 
For such relaxed boundary conditions, the existence of a variational principle requires to add 
an anomalous counterterm to the Einstein-Hilbert action. This counterterm arises as follows. 
Logarithmic divergent integrals at the spatial boundary of temporal hyperbolic initial and final 
slices in the action are regulated by introducing a hyperbolic radial slicing and by cancelling 
divergences at a finite radial hyperbolic cutoff. Thanks to the choice of hyperbolic cutoff, the 
regulation procedure preserves asymptotic Lorentz invariance. If this cutoff is pushed to infin- 
ity, the resulting regulated action will not be invariant under asymptotic radial diffeomorphisms 
since the action would be shifted by a finite piece. This finite action is precisely the anomaly 
and is proportional to the action for the first order boundary fields that one can infer from 
the bulk Einstein's equations. Therefore, the regulated action depends on the specific choice of 
hyperbolic radial foliation close to spatial infinity. A covariant phase space is defined only after 
one such slicing has been chosen and fixed. The anomalous counterterm is given by the anomaly 
multiplied by the logarithm of the radial hyperbolic cutoff. The action for asymptotically flat 
spacetimes without parity conditions has a radial hyperbolic slicing anomaly in the same sense 
that odd-dimensional anti-de Sitter spacetimes have a Weyl anomaly [21] 10 . The interpretation 
of this anomaly and the role that it may play in the quantization of asymptotically flat gravity 
largely remains to be understood. 

As a first step in the understanding of the structure of asymptotically flat spacetimes without 
parity conditions, we derived the asymptotic charges and the asymptotic symmetry group using 
the details of the Beig-Schmidt expansion. The hyperbolic radial slicing is not invariant under 
the action of (super/log)-translations. As expected from general considerations, the anomaly 
is however invariant under (super/log)-translations and Lorentz transformations. The relaxed 

10 The presence of anomalies is uncorrelated with the presence of central charges in the asymptotic symmetry 
group. Here, as in anti-de Sitter spacetimes in odd spacetime dimensions d > 5, the asymptotic symmetry group 
is finite-dimensional and not centrally-extended while anomalies are present. 



27 



boundary conditions lead to several additional non- vanishing conserved charges at spatial infinity 
such as logarithmic translation charges, supertranslation charges and boundary Noether charges. 
The logarithmic translation and supertranslation charges are canonically associated with bulk 
(super/log)-translations that are however not allowed to act on the phase space. The boundary 
Noether charges are associated with the boundary Lorentz symmetries of the anomaly. When 
parity conditions do not hold, the Lorentz charges are non-linear functionals of the asymptotic 
fields and therefore differ from the standard ADM and AD formulas [28, 31]. The standard ADM 
and AD formulas are restored when parity conditions hold. The non-linearities have been treated 
in detail in this work. Such asymptotic non-linearities appeared already for asymptotically anti- 
de Sitter spacetimes. In Einstein gravity, the charges are linear functionals of field perturbations 
around anti-de Sitter [31, 51, 52]. However, non-linearities may appear when matter fields are 
present, see e.g., [53, 54, 55]. 

It is interesting to remark that the presence of non- vanishing charges associated with super- 
translations in addition to Poincare transformations is also a feature of null infinity where super- 
translations along the null direction are also associated with non-trivial charges [56, 57, 58] 11 . 
For regular asymptotic fields, one expects that supertranslations charges should be conserved 
at infinite past times of future null infinity or at infinite late times of past null infinity where 
the news tensor vanishes. Indeed, at such late or early times the expression of [56, 57, 58] be- 
comes conserved and proportional to the first order electric part of the Weyl tensor and matches 
qualitatively with our expression (4.87). The shifts of Lorentz charges due to a change of radial 
hyperbolic cutoff are also reminiscient of the ambiguities in defining Lorentz charges at null 
infinity due to the action of supertranslations. It would be interesting to make that qualitative 
agreement more precise by comparing the precise definitions of supertranslations. 

We derived the asymptotic symmetry group using a generalization of the Poisson bracket 
taking into account boundary contributions to the charges. We obtained that the Poisson 
bracket is well-defined for all asymptotic symmetries: translations, logarithmic translations, 
supertranslations and Lorentz transformations. We derived the complete asymptotic symmetry 
algebra of conserved charges and we obtained that it is not centrally extended. The situation here 
can be contrasted to bulk infinitesimal diffeomorphisms which induce Killing symmetries and 
conformal Killing symmetries of asymptotically AdS spacetimes in odd dimensions as analyzed 
in [24, 25, 48, 62]. First, the dependence of the Lorentz charges, associated with Killing vectors, 
upon the choice of regulator is analogous to the shift of the stress-tensor by Weyl anomalous 
terms [24, 25]. In our case, logarithmic translations are also present and they are also shifted 
under a change of regulator. In anti-de Sitter, infinitesimal diffeomorphisms associated with 
boundary conformal Killing vectors are well-represented by the Poisson bracket even though 
they may act non-trivially on the action [48]. Indeed, the action only varies by a c-number 
which depends on the boundary conditions while the dynamical phase space is preserved. The 
non-conservation of the associated charges is related to this c-number. In asymptotically flat 
spacetimes, translations are also boundary conformal Killing vectors. Four-momenta as well as 
supertranslations are always exactly conserved and they do not vary under a change of regulator. 

No exact solution of vacuum Einstein's equations is known to us which breaks parity condi- 
tions. Such a solution would possess twelve boundary Noether charges in addition to Poincare, 
logarithmic translation and supertranslation charges. The boundary Noether charges are the 

llr The proposal of considering the so-called superrotations [59, 60] in addition to supertranslations at spatial 
infinity would also require a careful regularization of the super Lorentz charges, see [61], which would go beyond 
the considerations of this work. 
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Noether charges of the actions for the first order fields associated with the boundary Killing 
symmetries or equivalently with the asymptotic Lorentz Killing vectors. A subclass of those so- 
lutions exists as an analytic series expansion at spatial infinity. Indeed, one can consistently set 
the logarithmic terms in the expansions (3.14) and (3.25)-(3.26) to zero and still obey Einstein's 
equations by fixing six linear combinations of the boundary Noether charges to zero, see Section 
4.3 for details. Then, the original Beig-Schmidt expansion [35] which uses only polynomials in 
p is a consistent analytic asymptotic solution of Einstein's equations at all asymptotic orders 
which has six boundary Noether charges. We leave the existence, or not, of a regular solution 
in the bulk with such charges as an open question. 

As a side note, we could also generalize the construction of asymptotically flat spacetimes 
to include NUT charges. It would require to allow the field fc a j to contain Dirac-Misner string 
singularities. We present in Appendix D a preliminary lemma useful in that context. As a 
consequence of the lemma, the NUT four-momenta can be defined as the surface integral of 
the first order magnetic part of the Weyl tensor as (4.114), in complete analogy with the four- 
momenta (4.84), defined as an integral of the first order electric part of the Weyl tensor. The 
asymptotic symmetry group would be modified, see also [50, 63, 64]. Due to singularities in k a b, 
the definition of Lorentz charges would require a careful treatment of integrals on the sphere, 
especially the boundary terms on the sphere that we ignored in this work, which would make 
the analysis technically much more involved. 

We started this work by emphasizing the following puzzle: how to reconcile the two facts that 
parity conditions do not seem to be required for defining a variational problem for asymptoti- 
cally flat spacetimes while parity conditions are required for defining the covariant phase space 
symplectic structure and canonical two-form. We found that the puzzle can be resolved in the 
framework of regulating asymptotically flat spacetimes with hyperbolic cutoffs by considering 
the boundary terms in the variational principle on hyperbolic temporal cutoffs. When these 
boundary terms are considered, we showed that a variational principle which allows for asymp- 
totic Poincare diffeomorphisms exists only when parity conditions are imposed. This argument 
was argued to be independent on the particular counterterms chosen at spatial infinity, e.g., in 
the choice of the Mann-Marolf prescription. 

We stressed in the introduction that supertranslations are not on the same footing between 
the canonical formalism, where some supertranslations are allowed to act on the fields but are 
pure gauge, and the standard treatment of the covariant phase space where supertranslations 
are completely fixed. We partially resolved that tension by constructing a covariant phase space 
with parity conditions on the first order fields that admits some supertranslations that are pure 
gauge. This covariant phase space enlarges the one defined in [5, 2] by allowing to vary the 
dynamical field k a b (restricted to be traceless, divergence- free and regular). This phase space 
therefore admits non-trivial first order part of the magnetic Weyl tensor ^ while still 
having zero NUT charge. We showed that a sufficient condition to have a well-defined vari- 
ational principle on this phase space is to fix the trace and divergence of k a b while using the 
Mann-Marolf variational principle. Interestingly, supertranslations that leave the trace of k a b 
invariant are precisely the supertranslations that are conserved on the phase space. Also, these 
supertranslations are precisely the covariant analogue of temporal and radial supertranslations 
in the canonical formalim [1]. We only consider these supertranslations in this work. Logarith- 
mic translations and parity-odd supertranslations do not act on this phase space because of the 
parity conditions. Parity-even supertranslations act on the phase space but are pure gauge. This 
covariant phase space has all properties that one expects such as the Poincare group as asymp- 
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totic symmetry group. The anomaly, even though non-vanishing off-shell, does not contribute to 
the classical dynamics as a consequence of the parity conditions. Also, Poincare charges are lin- 
ear functionals of asymptotic fields and reduce to the standard expressions for conserved charges. 
An auxiliary question remains on the role of the so-called angular supertranslations (see Section 
3 for details), which can be summarized as follows. In the Hamiltonian boundary conditions 
stated in [1], angular supertranslations are allowed transformations and are pure gauge while 
these diffeomorphisms are not allowed in the covariant phase space where one fixes the mixed 
components g pa in the Beig-Schmidt expansion. The introduction of mixed components g pa in 
the Beig-Schmidt expansion would therefore be necessary in order to fully obtain the Lagrangian 
analogue of the Hamiltonian boundary conditions stated in [1]. We leave such a generalization 
for future work. 

The second puzzle that we raised in the introduction was related to the status of log- 
translations and parity-odd supertranslations: should they be considered as unphysical dif- 
feomorphisms (e.g. are they gauge transformations in an enlarged phase space) or physical 
diffeomorphisms (are they associated with conserved charges)? We attempted to resolve this 
puzzle by constructing a phase space where these transformations are allowed diffeomorphisms. 
In fine, we only constructed a collection of phase spaces without parity conditions where a radial 
hyperbolic foliation is required to uniquely define each phase space. Therefore, even for such 
generalized phase spaces, log-translations and parity-odd supertranslations are not allowed to 
act on the phase space because of the anomaly shift of the action under such a transformation. 
However, the diffeomorphisms are associated with non-vanishing conserved charges in general. 
The puzzle therefore remains. Lacking a clear interpretation of the covariant phase spaces with- 
out parity conditions, we leave to further work the elucidation of the role of these asymptotic 
diffeomorphisms. More generally, it is not clear what the anomaly tells us about the quantization 
of Einstein gravity. Nevertheless, since the anomaly can be derived from first principles from 
the Einstein-Hilbert action under reasonable assumptions, we think that it may find its place in 
(and maybe contribute to formulate) a future theory of quantum gravity in asymptotically flat 
spacetimes. 
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A Comparison of 3+1 and covariant boundary conditions 

The hyperbolic and cylindrical representation of spatial infinity are valid in the limits p — > oo 
and r — > oo, respectively. The key change of coordinates is the one mapping flat spacetime from 
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the hyperbolic to the cylindrical representation of spatial infinity 



t 2 t 
l-_ , T = arctanh(-) . (A.120) 



V r 

The hyperbolic and cylindrical representations coincide asymptotically in the limit where ADM 
time is kept finite, t/r — > which is equivalent to r — > 0. In that case, p ~ r asymptotically. 

In order to obtain the form of the metric in r, t coordinates, we expand the right-hand side 
of p, t in powers of t/r and we expand the Beig-Schmidt fields in Taylor series around r = 0, 

<T{r,e,<f>) = <r(eA) + - r ^(eA) + ^i(eA) + o{r^), (A.121) 
k ab (T,e,c/>) = k ab (e, ( f > )+ t -k: b (e, ( t > ) + ^ lab (e,^ + o(r- 3 ), (A.122) 

i ab (T,9,<P) = i ab (9, ( p) + t -il b (9^)+0(r~ 2 ), (A.123) 



r 



h^{rM) = h%(e,0)+ t -h: b i2 \e,0) + O(r- 2 ), (A.124) 

where we define a n (9, <j>) = <9 T a(0, 0, </>), = d T k ab {0,9,(f>), ~{{9,0) = d T d T a{0,9,(f>), 

lab(9A) = d T d T k ab (0,9,cl>), %l h = d T i ab (0,9,cl>), h^ b {2 \e,ct>) = d T h%>(0,9,<l>). We will keep the 
same notation for canonical fields in Hamiltonian formalism as fields in Lagrangian formalism 

a(0, 9, <j>) = a(9, 4>), k ab (0, 9, <f>) = k ab (9, 0), 
U0,M) = iab(0A), h%>(0,9,<t>) = hf b \9,<f>). (A.125) 

The tensors decompose into scalars, vectors and two-dimensional tensors under decomposition 
into temporal and spatial components. The meaning of the notation should be clear in either 
Hamiltonian or Lagrangian context. The fields j(9,(p) and j ab (9,<j)) are determined from the 
equations of motion of a and k ab . After a straightforward computation, we obtain 



1 2a a 2 + 2ta n . _ 2 , 
9rr = 1 + — + ~ 2 +o(r 2 ), 

(2) 

W = -~M-t " ^ ' S/ ^ rC ' T " +o(r~ 2 ), (A.126) 



- K T c t „ 

W» = r 2 5(t + (fc tC -2a 5( Jr + logr(i c J + (/ l ^+t^-2^ c J+ O (r ), 
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for the canonical fields and 



(detg^) 





{detg^c) 



(A.127) 



(detg (l ) 




for the conjugate fields. Here, we denote by (k — k terms) terms quadratic in k a b which do 
contribute to the finite part of the conserved Lorentz charges but that we omit here for simplicity. 

Let us finally discuss how the notions of parity are related between Beig-Schmidt fields and 
canonical fields. A field on the hyperboloid is parity-time reversal even if it is invariant under 
the combined transformation of inverting the hyperboloid time r —> — r and doing a parity 
transformation (6, (ft) — > (ir — , (ft + it) . Fields in canonical formalism are parity-time reversal 
even if their components in Cartesian coordinates do not transform under three-dimensional 
parity and if the components of their conjugate momentum in Cartesian coordinates transform 
with an overall sign under parity. From the dictionary of the Beig-Schmidt asymptotic fields 
in 3+1 decomposition, we see after switching from spherical to Cartesian coordinates that the 
even parity-time reversal conditions on a and k a b lead to parity-time reversal even first order 
canonical fields on the initial time slice t = 0. 

B Classification of symmetric and divergence-free tensors 

As discussed in the main text, conserved quantities associated with Poincare generators can 
be built out of symmetric divergence-free and traceless (SDT) tensors, or more generally out 
of symmetric and divergence-free (SD) tensors. It is useful in order to establish unicity of the 
definition of conserved charges and in order to understand the structure of the linearization 
stability constraints at second order to provide the classification of all possible SD tensors that 
one can built up from quadratic terms in a and k a b and their derivatives. In fact, for both 
purposes, it is sufficient to classify SD tensors whose curl are non-zero since the integral on the 
sphere of a SD tensor with zero curl contracted with a Killing or conformal Killing vector is 
identically zero (see Appendix B of [12]). 

All SD tensors built out of quadratic terms in the first order fields can be formed from 
symmetric tensors M a b obeying T> b M a b = T> a M and which we call tensor potentials. A complete 
set of SD tensors consists of SD tensors given by K a b = M a b — Mh^j and of SDT tensors obtained 
by acting with successive curls on M a h or equivalently by acting with successive symmetrized 
curls on K a b. Indeed, an SDT tensor can be constructed from T a b = e^VcM^b = ^ c d(aP c K b) • 
As we have just emphasized, the curl of a tensor potential might be trivially zero. Therefore the 
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SD and SDT tensors whose curl are non-zero will turn out to be classified using the equivalence 
of classes of tensor potentials where two tensor potentials are equivalent if their difference has 
a trivial curl. We will denote one representative of such equivalence class of non-trivial tensor 
potentials a RNT tensor potential. 

We explain in the next subsection the procedure we have followed in order to obtain and 
prove that we have obtained all RNT tensor potentials needed to construct any SD or SDT 
tensor whose curl is non-zero. 

B.l Algorithm for classification 

The first order fields a and k a b obey decoupled linear equations. We thus consider separately the 
quadratic combinations (a, a), (k, k) and (a, k). For each case, we use the following procedure: 

1. We start by listing a basis of rank two symmetric tensors with m derivatives built out 
of quadratic terms which are independent on-shell. It exists a number of derivatives m* 
such that for all m > to* the number of terms in that basis is maximal, i.e. no new 
tensor structure appears at higher order. At lower values m < m* not all possible tensor 
structures can appear due to a lack of derivatives. We find to* = 3 for (a, a), to* = 3 
for (k, k) and m* = 4 for (k, a) tensors. Due to the presence or not of the epsilon tensor 
depending if m is even or odd, the general form of a rank two symmmetric tensor built 
out of linear combinations of the basis takes a different form. We denote this tensor as 
Q^ab^ or Q^ab +1 ^ an( ^ we P roy ide its general form. 

2. We continue by deriving a bound on the possible SDT tensors that one can build at a fixed 
number m > to* of derivatives. We simply compute the number H of linearly independent 
tensors which obey both T^ b Q^ = and Q^ " = where equalities here are valid 
up to terms with lower derivatives. At this stage, the number H is only a bound on the 
number of SDT tensors at order m because none of them has been constructed fully yet. 
We obtain that H = 1 in the (a, a) case, H = 3 in the (k, k) case and H = 2 in the (a, k) 
case for both m even or odd. 

3. We then derive the explicit form of all RNT potentials, SD tensors and SDT tensors at 
each low value m < m* of derivatives by enumeration. We write a basis of rank two 
symmetric tensors built out of quadratic terms which are independent on-shell with at 
most m derivatives for each m < to* and impose the RNT or SDT conditions. The 
SD tensors K a b, whose curls are non-zero, are obtained from the RNT potentials by the 
correspondence n a b = M a b — h^M^. 

4. We finally observe that there are exactly H SDT tensors which have at most m* derivatives 
and at least one term with to* derivatives. This provides a proof that each candidate SDT 
tensor exists at order to*. We then note that the SDT tensors obtained by acting with 
the curl operator on these tensors form a basis for SDT tensors at order m* + 1 and by 
successive iterations at each order m> m*. Since there are H SDT tensors at each order 
m > to*, there cannot be any other SD tensor which is not traceless but whose curls 
are non-zero or equivalently any RNT tensor at order to. Otherwise, there would be one 
additional SDT tensor at order m + 1 by applying the curl operator but this would raise 
the number of SDT at level m + 1 to H + 1, which is not the case. 
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5. We conclude that all RNT potentials and SD tensors whose curl are non-zero are classified 
by the explicit tensors that we build out of terms with up to m* derivatives. At higher 
order m > m* in derivatives, all SD tensors whose curl is non-zero are traceless and can 
be obtained by applying curls on the RNT potentials. 



B.2 (a,a) SD tensors 

This analysis was already performed in [12]. Let us summarize the results here in our new 
notation. We have m* = 3 because the tensor structure t a ObcPd first appears at order m = 3 
and no other tensor structure appears at higher m. There are no tensor potentials at order 
m = and m = 1. At order m = 2, there are two independent tensor potentials 

M^ I] = (5a 2 + a c a c )h^ + 4aa ab , 

M% a ' a ' II] = (V a V b + h^)a 2 . (B.128) 

The first one is a RNT potential whose successive curls generate the unique SDT tensor at each 
order m > m*, while the second one has a zero curl. The two SD tensors associated with the 
tensor potentials are given by 

M [2 *,*,!) _ M [2,«,*,r] h (o) 

ab ab 

(2a 2 - 2a c a c )h a ° b ] + Aaa ab , (B.129) 

M [2 a,a,II] _ M [2,a,a,II] h (0) 
ab ab 

2a a a b + 2aa ab + h$ (4a 2 - 2a c a c ) . (B.130) 
The second SD tensor has a zero symmetrized curl. 



v ab 



K 



[(7,(7,11] 

ab 



B.3 (k,k) SD tensors 

The classification of (k,k) structures is considerably more involved than the classification of 
(a, a) structures. In order to simplify the identification of a basis of independent tensors on- 
shell, we will make an efficient use of the relations 

V [a B^ c = , V [a k b]c = -e abd d . (B.131) 

We start by listing the independent structures quadratic in and its derivatives. Then, we 
add an independent subset of structures of the form (B^\ k) such that no linear combinations 
are of the form (B^\ B^) and eventually we add a subset of independent (k,k) structures 
such that no linear combinations are of the form (B^\ B^) or (B^\ k). To look if such linear 
combinations exist, we just need to take into account the equations (B.131). 
For an odd number (2n + 1) of derivatives we find the general form 

qT l) = 

a e^V^.-V^V^B^ ce V il ...V in _ 1 B^ d + b e cd(a V h ...V^V^V^V^ ...V^V f k d e 
+eV i \..V i ^V c V d B a l b ) V il ...V in _ l k cd + fV i K..V i ^B^ ) V h ...V in _ 1 V c V {a k b f , 
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while for an even number (2n) of derivatives we find 

aV h ...V in _ 2 V c B d 1 Jv h ...V i ^V c B^ de h { $ + b V il ...V in _ 2 V c B^V i \..V i ^V c B^ d 
+c V h ...V in _ 2 V c V d B ( ^V i \..V i ^B^ cd + d e^V^.-.V^Vf V c k b p h ...V in _ 2 V f B^ d 
+e e a n a V il ...V i »-*VfV b) BW c V h ...V in _ 2 V f k de + f V h ...V ln _ 2 V e V {a k cd V^ ...V^V e V b) k cd 
+gV il ...V in _ 2 V e V {a V b) k cd V i \..V i ^V e k cd + hV H ...V tn _ 2 V c V d k ef V i \..V^V c V d k^ . 

(2) 

We deduce that to* = 3. Indeed, when to = 2 the third term in Q ab does not exist while 

when m = 3,4 or higher all terms in Q^P exist. Looking at m > to* and imposing the SDT 
condition, we find after a straightforward analysis that there can be at most three independent 
SDT tensors. We thus have H = 3. 

We now need to construct RNT potentials and SDT tensors at each order m < to*. At 
m = 0, there are no tensor potentials and no SDT tensors. At m = 1, we have 

Off = « M.^, + & ^(% C + c BS ) fc a, fc2 ) > (B.132) 
and one easily checks that there are no SDT tensors but there is one RNT potential : 

M [ ^ ] =-4B^ + B^h^. (B.133) 

At m = 2, we have: 

Q$ = « B$BV h <® + b bVb§> c + c BW rf + d e cd{a V b) B^ '** 

+e V {a k cd V b) k cd + f k cd V {a V b) k cd + g V c k de V c k de h^ 

+hk c {a k b)c + ik cd k cd h^ , (B.134) 

where we also introduced the structures with to = derivatives. Here, we get: 

Q( 2 ) = [3a + 6-2c] J B( 1 ) c ^ ) +P 6 A; cd P 6 A: cd [ e + 3 5 ] + [3/ + / i + 3^ cd fe cd , 
P a Q( 2 ) = [6a + 26 - 4c]BW cd V a B^ + [2e + 6 5 ]P a P fe fc cd P b fe cd 
+[6f + 2h + 6i]k cd V a k cd , 

V b Q<$ = [2a + b-2d}B^ cd V a B^ + [^-^}e cda V% e V b B^ d 

+ [2c - 2d - 4e + 4/ + 2/i] e cda k ce B^ d + [ e + f + 2g]V a V b k cd D b k cd 
+k cd V a k cd [e + 7f + h + 2i] , (B.135) 

where we made use of the relations: 

^□P a fc cd = 5k cd V a k cd + Ak cd V c k ad , k cd V b V c V d k ab = 7k cd V c k ad , 
k cd V b V a V b k cd = k cd UV a k cd - 2k cd V c k ad , k cd V c k ad = 2e cda k\B { p d + k cd V a k cd . 

We obtain that tensors Q a 2 b satisfying V b Q a 2 b = V a Q^ are of the form 

rrnY^ + m 2 M [2 ' k ' k ' I] + m 3 Mj Wil + m 4 M [2 ' k ' k ' IH] , (B.136) 
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where 

= -4 5«Bg> c - 2 e cd(a V%pV * _ 2 ecd{a V b) B^ c k de = e cd(a V c M^ d 

M[f AI] = l -k cd k cd h { ^ - k ac k% + l -V c k de V c k de h^ - l -V a k cd V b k cd , 

M [2k,k,II] = A B W B Wcd h (0) +B w B wc 

ab a ca ab c(a b) ' 



M [2,k,k,ni] = l _ k cd V{J)b)kcd + \ V{ak - d V h) k cd + ^k cd k cd h { 2 . (B.137) 



We thus see that is the unique SDT tensor obtained from the RNT potential M a \' k ' k ^ , that 

M a 2 b ' k ' k ' I] and AfjJ*'*'" 1 are two new RNT potentials and that M [2 b ' k ' k ' HI] is a tensor potential 

which has vanishing curl as it is of the form {V a V b +h a ° b )k cd k cd . From the three tensor potentials 
found, we can define 3 SD tensors: 

= lk cd k cd h a ° b ] - k ac k c b + ^V c k de V c k de h^ - ^V a k cd V b k cd , 

n [Kk,III] = ^ M %k,k,III\ _ M [2,k,k,III] h {0)^ 

= (-^V c k de V c k de - h^h^ + ^P (a & cd P b) fc cci + ^k cd V (a V b) k cd .(B.138) 
Now, at m = 3, we obtain: 

q£> = ae cd( J) fe) 5( 1 ) ce i?( 1 ) d + &e cd ^^ 

+d V c B d X }v c k d ^ + e V c V d B a X b \^ + / B^V c V {a k b d , (B.139) 

up to terms with lower derivatives. One can explicitly construct three independent SDT ten- 
sors which can also be obtained as curls of the three previous RNT potentials. We have thus 
completed our algorithm. We have three towers of SDT tensors generated by the three RNT 
potentials M% k ' k] (which leads to y} 2) ), M [2 b ' k ' kJ] and M [2 b ' k ' k ' H] . 

B.4 (a, k) SD tensors 

For the (a, k) case, we find that a generic tensor with 2n + 1 derivatives is of the form 

+c e cd{a V i K..V i -- 1 V b) k%a\^ in _ 1 + d e cd{a V*K..V^k c *a b d ^...^ 
+e V il ...V in _ 1 B^a cdi ^ i ^h a ° b ) + / ^...D^^*"^ , (B.140) 
while, for an even number 2n of derivatives, it is of the form 

QaT ] = aa cdh - i ^V h ...V in _ 2 V c V {a k b)d + b a cdil - in - 2 V h . . .V in _ 2 V {a V b) k cd 
+ c a cdi ^ i ^V il ...V in _ 2 V c V d k ab + d a cdh _.. in _ 2(a V i \..V i ^V b) k cd 
+ e a cdll .. An _ 2(a m...V i ^V c k b d + / a abcdn ... ln _ 2 V l \..V i ^k cd 

+ gh ( ^ (T cdil ... in _ 1 V i \..V i ^k cd . (B.141) 
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(3) 

We find that m* = 4. Indeed, for m = 3 derivatives, the last term in Q ab does not exist while 
for m = 4, 5, . . . all terms in Q a ™^ exist. For any m > 4, one can check that there are at most 2 
SDT tensors. We therefore find H = 2. At lower levels m = or m = 1, we see that there are 
no SDT tensors and no tensor potentials. At m = 2, we have 

Qab = a ak ab + b e cd(a k d b) a c + caB^ + da c V c k ab 

+e a c V [a k h)c + / a c(a k b) c + g a cd k cd h a ° b ] . (B.142) 

There is no SDT tensor, but there is one RNT potential 

M a 2 b aM = °Kb - a c V c k ab + a c V (a k b)c + a c(a k b f - ^ a^h^ . (B.143) 
Its curl gives an SDT tensor 



l<, 'n.k'i, - - a' V. lC, ■ (B. H I) 



At m = 3, the general SDT tensor or RNT potential can be written as a linear combination of 
a basis of terms with 1 and 3 derivatives (terms with 2 derivatives decouple) 

Qab = « Sg% c + b hatred + c M^f*! + * e cd(a V b) k d e a ce + e a e V e B^ 

+fe cd{a k%a c + g oB% . (B.145) 

We find one SDT tensor which is obviously Z a 3 b and a new RNT potential M a l' a ' k] 



As expected, at m = 4, one can check that we have 2 SDT tensors. The algorithm is therefore 
completed. The two RNT potentials that generate the two independent towers of SDT tensors 
are M^ h ' a ' k ^ and M n 3 h ' a ' k ^ . To each of these potentials corresponds a unique SD tensor 



^ab lv± ab "'ab 



o,k,I] = M [2,*,k] _ h (0) M {2,a,k] 

1 

2 

>,< 

( 3 ) iOrW^ C jinW 1,(0) 



(rk ab + a c(a k b) c - la cd k cd h a b ] - a c V c k ab + a c V (a k b)c , (B.147) 



K [a,k,II] = M ^M_ h m M [3,a,k] 



= Z^ + 2B^a b f-a^B^h^ + 2aB^. (B.148) 

All other SDT tensors, with m > m*, are then constructed from linear combinations of successive 
curls of those SD tensors. This ends the classification of RNT and SDT (a, k) tensors. 
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B.5 Properties 

In order to simplify the linearization stability constraints given in Section 4.3, we would like to 
understand what currents K ab associated with SD tensors n ab constructed in the previous 
subsections can be expressed as total divergences. Here, we show that the currents associated 
with the two independent (a, k) SD tensors and two linear combinations of the three (k, k) 
SD tensors which are not traceless are total divergences. We are thus left with two currents, 
one quadratic in a and one in k, which can be identified as the currents constructed using the 
stress-energy tensors of the counter-term actions and up to total divergences. 
For the (a, k) SD tensors, one can first check that the current: 

K hM^ =v b^_ ^ ^ + v%[a ah]c + e aV[ah]c + e (R149) 
can be expressed as a total divergence and will thus not contribute to the charges. As we know 

(3) 

from [12], this also implies that its symmetrized curl, the SDT tensor 2TV, will not contribute 

to the charges. If we look at the definition of given in (B.148), we are left with checking 

that 

K(% C - o^BS fcff + laB^y = 2V b {D^ [a aB® - t [a B%K e + f ^fig) (B.150) 

to see that k^' 1 will not contribute either to the charges. This proves that no SDT or SD 
tensors quadratic in (a, k) will contribute to the charges. For the terms quadratic in (k, k), one 
can show by inspection that: 

V h (2i c k d[a V%f - V c i d k c [a k b] d ) 

= (v c k d(a V d k c b) + Ak c(a k b f - k cd V c V d k ab Y 

= ( 7 bWbU h <fi) _ 6 flWsg) ■ + 4 ecd(aV% e B (l) 4 + 4 €cd{aVb)B p c kd e 

+ V (a k cd V b) k cd - k cd V {a V b) k cd + Ak c (a k b)c + k cd k cd h a Q b ] Y 
= ( - 2^i? - ~ 14$F* - , (B.151) 

and also: 

v b (v c i [a k b d k c d - e [a v c k b d k cd - e k\ a v h] k cd + e k cd v [a k b d ) 

= ( - 2 k\ a k b)c - \ k^h® - \v c k de v d n[ Q b ] 

+V c k d a V b) k cd - V c k d(a V c k d } + k cd v c v {a k b)d Y 
= ( - 5 B% B V <* h<$ + 6 B«B« c _ 2 ecd(aV% e B W d _ 2 ecd{aVb)B W c kde 

+ k cd V {a V b) k cd - l - V c k de V c k de h a ° b ] - 2 k\ a k b)c - I k^h^Y 
= (Y<Z> + 2^ + 10«S^ + ^ III] Y • (B.152) 
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The SD tensor k^' 111 ^ has a trivial symmetrized curl. It is thus associated with trivial charges 
when contracted with a Killing vector. This result is also reflected in the fact that the current 

4f "V = -\v b (2k c % a V b] k cd + h cd k cd V [a £ b] ) (B.153) 

can be written as a total divergence. Eventually, we notice from (B.151) and (B.152), that the 
current constructed out of K^ k ' can be written as a total divergence and that the equality 

Y m^)b = _ 2K [k,k,i]^ )b (B154) 

is true up to a total divergence. We have thus shown in particular that the charges constructed 
from the current n ab £? q n where n ab is defined in (4.49) can be written as 

<j> d 2 Sy^h(0)K ab £f 0) n h = I d 2 Sv^hW (47>'] + l -Y^) e {0) n h . (B.155) 



C Equations of motion in Beig- Schmidt form 

In this appendix, we will first review the 3+1 split of Einstein's equations G pv = for the four- 
dimensional metric g pu on three-dimensional hypersurfaces h ab which have a spacelike normal n^, 
n^n p = +1. We will then expand those equations asymptotically using Beig-Schmidt coordinates 
up to second order in the expansion while keeping a ^ 0, k ab / and i ab ^ 0. This extends all 
previous treatments, see e.g. [35, 36, 16, 10, 3, 13, 14]. We will eventually rewrite the second 
order equations in terms of two, much more compact, equivalent systems. 

Our notation is as follows. We use greek letters for four-dimensional quantities and 

latin letters a, b, ... for three-dimensional quantities. The letter p always denotes the radial 
coordinate that we normalize as n^d^ = (1 + ^) 1 -§^- We define the metric h pu = g pu — n p n u . 
The extrinsic curvature is defined by K ab = hah^V^n^ and its trace is K = h ab K ab . In our 
case, it reduces to K ab = ^n p d p h ab . 

Covariant derivatives associated with g pv , h ab , and are respectively denoted by V M , D a 
and T> a . We have D a = h a ^V „. The four-dimensional curvature, Ricci tensor and scalar are 
denoted by R% p(T , R pu and R while their 3d counterparts are denoted by 7Z a bcd ,1Z ab and TZ. 
We use the Misner-Thorne- Wheeler conventions R a bcd = R c f b = d c F bd + T a C( £ bd — (cod), 
Rab = R c acb - 



C.l The 3+1 split 

Einstein's equations can be split into a set of three equations when appropriately projected along 
perpendicular and/or normal directions to the hypersurface. This provides us with Hamiltonian 
and momentum equations of motion (these equations contain time derivatives and therefore are 
not constraints) and equations of motion on the 3-dimensional hypersurface which respectively 
read 

H = -2G pu n p n u = K + K ab K ah - K 2 = , 

F a = h a tt n u R fai = D b K b a — D a K = , (C.156) 
F ab = KK R w = K ab - (C n K) ab + D {a a b) - a a a b - KK ab + 2K a c K cb = , 
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where C n stands for the Lie derivative along the unit normal and a p = n p V p n p is the 4- 
acceleration. 

The Hamiltonian equation can be further simplified by taking the trace of the equation of 
motion so that we get 

H = -C n K - K ab K ab + D a a a - a a a a = . (C.157) 

Using our values for the shift N a = 0, we get D^ a a b ^ — a a a b = —N~ l D a D b N and the equations 
reduce to the ones of Beig and Schmidt [35] 

H = -C n K-K ab K ab -N- l h ab D a D b N = 0, 
F a = D b K b a — D a K = , 

F ab = n ab -N- 1 d p K ab -N- 1 D a D b N-KK ab + 2K a c K cb = 0. (C.158) 
C.2 Equations of motion in the radial expansion 

We now expand these equations using our boundary conditions (3.14). The inverse metric is 
expanded as 

h ab (p, x c ) = p- 2 h^ ab - p~ 3 h^ ab - \npp- 4 i ab - p~\h^ ab - /i«>« c6 ) + 0(p~ 5 ) . 
The extrinsic curvature admits the simple expansion 

K ab = ph$ + ( h% - <rh<$) + \ ( U ab - \h[% + o> h <$) + 0( P - 2 ) , (C.159) 



and we also have 



K<b = ^-^-^>+^(-^^ 



-In 



The covariant derivative requires an expansion of the Christoffel symbols 

rfc = r (0) a bc + P - l r^ \ c + \n PP ^T^l a + a bc + o( P - 3 ) , (c.ieo) 

where 



v {ln,2)a = ^^a^^a^^^^ (C161) 



be = j V C b + b c ~ bc J ~ 2 \ Ucn db f U V l dc - u d% c f 

The expansion of the three-dimensional Ricci curvature tensor is 

n ab = k[v + p-^2 + inpp- 2 n^ + p -2 U m + 0(p - 3) (c 162) 

The zeroth order Ricci tensor is the one constructed with the metric h„°? . The first order Ricci 
tensor and the tensor 1Z are 
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n 



(i) 

ab 



,(l)c 
ab 



-v h 



r (l) 



P c P 6 fr& + V c V a h£) - V a V h hV - V c V c h 



n 



(ln,2) _ 1 



P c P b i ac + P c P a i 6c - P a P b i - V c V c i ab 



ab 2 

and the second order Ricci tensor reads as 

n 



,(2) _ 1 
ab 2 



V°V b hW + P C P«4? - V a V b h™ - v c v c h { 2 



-\v d 



hM*(VJ$ +V h h$ -V c h 



(1) 
ab 



Finally, the equations can be expanded as 



H = p~ 3 H<U+ In pp- 4 H^ + p- 4 H^ + 0(p- 5 ), 
F a = P- 2 F^ + In PP -*F^ + p-*FP + 0(p-*) , 



(C.163) 



ab 



At zeroth order, we only have F 



(o) 



ab 



K 



(o) 

ab 



2h 



(o) 

aft 



which implies that the boundary metric 



is three-dimensional de Sitter spacetime. At first order, the Hamiltonian equation = is 
simply 



(□ + 3)cr = 



The momentum equation F^ = is 



V b k ab = 



and the radial equation of motion F^ is 



(□ - 3)k ab = 



(C.164) 



(C.165) 



(C.166) 



after we set k = 0. At second order we easily get for the logarithmic terms i?(' n > 2 ) = 0, Fa 



(Zn,2) 



and F, 



(ln,2) 



ab 



o, 



V b i ab = 0, (□ - 2)i ab = 



(C.167) 



For the finite terms at second order we find 

2 4 



+9a 2 + aV 2 a + h^ ab V a V b a + V b aV a h ( a \ ) - ^V a aV a h^ . (C.168) 
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Using only a and k ab = h^ b + 2ah® b , and also k = i = 0, we obtain 



h® = 12a 2 + a c a c + \k cd k cd + k cd a 



cd 



where we also made use of the first order equations of motion. We also have: 

= V b K^ b - V a K^ + T^K^ - F^ c Kp b = , 
which amounts, after simplifications, to 

Vbh< ab = \V b k ac k b c + V a (a c a c + 8a 2 - \k cd k cd + k cd a cd ) 



(C.169) 



(C.170) 



(C.171) 



The radial equation of motion can be obtained after a straightforward, although tedious, com- 
putation and we find the quite intricate form 



(□ - 2)b$ = 2i ab + NL a6 (a, a) + NL o6 (a, k) + NL afc (A;, k) 



(C.172) 



where the non-linear terms are given by 

NL o6 ((T, a) = V a V b (5a 2 + a c a c ) + h$ (-18a 2 + 4a c a c ) + 4aa ab , 

NL ab (a, k) = V a V b (k cd a cd ) - 2k cd a cd h^ + Aak ab + 4a c (D {a k b)c - D c k ab ) + 4a c(a k' 

NL ab (k,k) = k ac k c b + k cd (-V d V {a k b)c + V c V d k ab ) 



^D b k cd V a k cd + V a k c(a V b) k d c + V c k ad V c k a b - V c k ad V a k% . 



cud 



(C.173) 



Using the relation a c a abc = a c a cab + a a a b — h a ° b a c a c , one can rewrite the NL ab (a, a) non-linear 
terms as 



NL a6 (<r, a) = 6a c a c h a ° b ] + 8a a a b + Uaa ab - 18a 2 h a ° b ] + 2a ac a% + 2a abc a c 
The equations of motion reproduce the expressions of [35] when k ab = i ab = 0. 



(C.174) 



C.3 Equivalent systems for second order equations 

Here, we would like to rewrite the second order equations in terms of two tensors V ab and W ab 



V, 



— ~h ab + ~ji a b + Qab > 



IT,,, ee e*D c (h$-~i db + Q. 



w 

db I ! 



(C.175) 
(C.176) 



where are appropriate quadratic terms in (a, a), (a, k a b) or (k a b, k a b) that we will construct 

herebelow. We require that the V, W are SDT tensors that obey the following duality properties 

W ab + e a cd D c V db = K% , V ab - e a cd D c W db = -2i ab + K v ab , (C. 177) 
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where K^ b W are non-linear terms quadratic in a and k ab which are also SDT. Applying the curl 
operator on both equations (C.177) we obtain that V ab and W ab obey 

(□ - 2)V ab = -2i ab + K\ + e a cd D c Kj b , (C.178) 
(□ - 2)W ab = -2j ab + K% - e a cd D c Kl, (C.179) 

where j ab = —curl(i) ab . Our construction of the non-linear tensors Q ' , K ' goes as follows. 
In order for W ab to be traceless, we require to be symmetric. Using the Hamiltonian and 
momentum equation of motion, one can rewrite the symmetry condition of W ab as the following 
equation on Q\\ 

V b Q% - V a Qf b = -h bc V b k ac + V a (V + h cd k cd ^j . (C.180) 
The divergence-free conditions of W ab can then be rewritten as 



e a cd V c (v b Qf b + \v e k d} kA = , (C.181) 



which is a consequence of the previous equation. The equation (C.180) can be solved up to the 
ambiguity of adding to Q^ b tensors obeying T> b M ab = T> a M. We will fix the ambiguity in defining 
Q^ b since we would like to find one equivalent formulation of the equations of motion, not all 

f ab 

obtain 



possible formulations. By choosing a Q\X with the smallest possible number of derivatives, we 



QZ = ( -2a 2 + Y Q k cd kA - \k ac k\ . (C.182) 



Using again the Hamiltonian and momentum equation of motion, one can rewrite the traceless 

?o6 



and divergence-free conditions of V ab as the following equations on QY 



Q Va = \2a 2 + a c a c + - A k cd k cd + k cd a cd , (C.183) 

V h Qab = \v b k ac k b c + V a (a c a c + 8a 2 - h cd k cd + k cd a cd ^j . (C.184) 

This system has a unique solution up to the ambiguity of adding an SDT tensor to Q v b - We 
will do a definite choice for the ambiguity in defining Q v b as well. In fact, the SDT tensor K^ b 
can be computed using (C.177) and the equations of motion (C.173) as 

K v ab = -NL ab (a,a)-NL ab (a,k)-NL ab (k,k) + Q v b -(n-3)Q a v b 

+V a V%h$ + OZ) ~ h a ° b \hW + QY a ) (C.185) 

We will choose to fix the ambiguity of adding SDT tensors to Q v b by requiring 

Kl = 0. (C.186) 

After a tedious computation, we obtain simply 

Q v ab = {Qa 2 + a c a c + \k cd k cd + \v c k de V c k de )h^ + 2aa ab -2a a a b 

o o 

+4 ak ab - 4 a c V c k ab + 4 a c V {a k b)c + ia c{a k b) c - a^h^ 

-\KcK - *-V a k cd V b k cd + h cd V {a V b) k cd + , (C.187) 
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where Y\ is an SDT tensor given in (B.137). Remark that to perform this computation, one 
can separate the analysis of non-linear terms for each set of quadratic terms (k,k), (a,k) or 
(k, k) independently since those terms never mix in the equations. 
Using then the definition of in (C.177) we find 

K\\ = curl(M) ab = e a cd V c M db , (C.188) 

where 

M ab ^QZ + Q V ab (C.189) 

is a tensor obeying T> b M ab = T> a M. Using the classification of such tensors in Appendix B, we 
have explicitly 

M ab = M^ I] - M^ II] + 4m£ ° M + + M^ kM + M^ k ' IH] , (C.190) 



where 



M^ I] = (5a 2 + * c a c )hW + Aaa ab , 

M%° M = ak ab - a c V c k ab + a c V [a k b)c + a c(a k b) c - i a cd k cd h { ° b \ 

M a 2 b h ' k ' I] = ^rt^-A^ 

[2,k,k,III] _ t>. .LhQAucdi 



= i-Av a v b +h^)k cd k cd . (c.191) 



In summary, the equations of motion can be written in the form: 

w a a = v b w ab = o, 

(D-2)W ab = curl(2i + M) ab , (C.192) 
% a a = V b i ab = 0, 
(D-2)i ab = 0. (C.193) 

Using the curl operator, the first set of equations lead to 

Va = V b V ab = 0, 

(D-2)V ab = curl(curl(2i + M)) ab . (C.194) 

Note that in the main text we have preferred to use the symmetrized curl of n ab = M ab — Mh a ° b . 
This is indeed equivalent to the curl of the tensor potential M ab . One can also get rid of any 
SD tensor whose symmetrized curl is zero as a consequence of the properties of integration on 
the hyperboloid. 

D Lemma on singular tensors 

This lemma was derived in collaboration with A. Virmani. 



44 



Lemma 1 On the three dimensional hyperboloid, any scalar <I> satisfying □<& + 3<I> = defines 
a symmetric, traceless, curl-free and divergence-free tensor T ab = D a D^ + h^?^ which can be 
written as 

T ab = e a cd D c P dh) (D.195) 



where P ab is a symmetric, traceless tensor of the form 



Ai=0 



(D.196) 



where P^ 9 is regular and k a ^ are four singular tensors listed here below. 

The regular tensors P^ 9 can be deduced from a Lemma stated in [35] and referred to as Lemma 
2 in [12]. The four singular tensors k}$ can be derived by integrating equation (D.195) for 
$ = where are the four solutions of (□ + 3)$ = which are even under parity-time 
reversal and contain the harmonics I = or I = 1 on the two-sphere. Explicitly, 



C(0) 
C(2) 



cosh 2r 
coshr 



C(i) 



„ . , tanhr 

2 sinn r H | cos I 

coshr 



tanhr , 

2 sinn r H — sin cos <p Q3) 

coshr ' 



2 sinh r + 



tanh r 
coshr 



The four singular tensors can be written in the traceless gauge hf^k 



(or(M)ab 



(D.197) 
sini9sinc/>(D.198) 
as 



k(l)ab 



z {2)ab 



z (3)ab 













2 k— cos 6 
cosh r_ 
g ^ nn r cos26>-4fccos6>+3 



\ o fc— cos 9 c^tiV, ^ cos 29— Ak cos 0+3 

\ Z coshr 2 sin6» 

/ 



o tanhr - 2 a -8fc+9 cos fl-cos 30 „ n „t, _ 

\ ^coshr 3111 17 4 sine cosn r 





2 sine 







sin 2 i 

cosh r 



-8fc+9cose-cos3e t 

4 sine cosnr 



V 3 



V3 



tanh j 
cosh 7 



S^sincA 
cos # sin # cos i 



8fc _9cose+cos3e cogh r gin ( 



(D.199) 



3^ cos 5 sin& cos, 



i tanhr 
cosh r 



COS < 



tanh r 
cosh r 



cos t> sin W sin < 



2 sin J 6 

cos 4 e-4fccose+3 h ^ 

"3^7 cos 
-8fc+9cos6>-cos3g cogh COS 

2 sni J e 

cos 4 e-4fccose+3 cosh T gin - 



-4fecosg+3 cos hr COS. 



sir?" 9 

■, +9c oBg-coB30 coshr sin, 



2 sine 

3^cos#sin#sin</> 

cos 4 e-4fc cose+3 cosh sin , 

sin^ e 
8fc-9cos0+cos 36 
2 sin e 



cosh r COS i 



These tensors are regular in the north patch upon choosing k = +1 and in the south patch upon 
choosing fe = — 1. They are tranverse and obey the equation 



(□-3)*; 



(fi)ab 







(D.200) 



outside of singularities. The singular transition function between the south and north patches 
can be written as 
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&k(i) ab = k^) ab \south 



&k(2)ab = k(2)ab\south 



<%3)a6 = k( 3 ) ab \south 



These transition functions obey 

D [a 5k Mb]c = 0, (□ - 3)5k Mab = 0, h^ ab 5k Mab = 0, D b 5k Wab = 0, (D.201) 

on the hyperboloid outside the singular region 9 = and 9 = ir and obey the normalized 
orthogonality relations 

/ Sk^D^ = -tor 5 {li){v) , = 0,... 3, (D.202) 

Jo 

where C(^) are the four solutions of V a V b (^ + /i^C(» = given by 

Qo) = — sinhr, £(i) = cosh r cos 0, (D.203) 
C(2) = cosh r sin cos C(3) = cosh T sin & 4> (D.204) 

which are odd under parity-time reversal and normalized such that Co)<9 p + P~ 1 d a (^d a + 
o(p^ 1 ) = dfj, where <9 M = d t , di. 
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